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Abstract. We consider a non-trapping n-dimensional Lorentzian man- 
ifold endowed with an end structure modeled on the radial compactifi- 
cation of Minkowski space. We find a full asymptotic expansion for tem- 
pered forward solutions of the wave equation in all asymptotic regimes. 
The rates of decay seen in the asymptotic expansion are related to 
the resonances of a natural asymptotically hyperbolic problem on the 
"northern cap" of the compactification. For small perturbations of 
Minkowski space that fit into our framework, we show a rate of decay 
that improves on the Klainerman-Sobolev estimates. 



1. Introduction 

We consider an n-dimensional Lorentzian manifold endowed with the 
end structure of a "scattering manifold" motivated by the analogous def- 
inition for Riemannian manifolds given by Melrose [15J. Our manifolds 
come equipped with compactifications to smooth manifolds-with-boundary, 
i.e., we will consider the Lorentzian manifold (M°,g) where M is a m.w.b. 
with boundary X = dM. The key example is the radial compactification of 
Minkowski space , where X is a "sphere at infinity" with boundary 

defining function p = (\x\ 2 + t 2 + 1) _1//2 . 

On M the forward and backward light cones emanating from any point 
q £ M° terminate at dM in manifolds S± independent of the choice of 
q; we call S± the future and past light cones at infinity, and they bound 
submanifolds (which are open subsets) C± of X = dM, which we call future 
(C+) and past infinity (C_). In the case of Minkowski space C+ and C_ are 
the "north" and "south" polar regions (or caps) on the sphere at infinity. 
Further, there is an intermediate ("equatorial," on the sphere at infinity 
in the case of Minkowski space) region Cq which has as its two boundary 
hypersurfaces 5+ and S-. We assume that the metric g is non-trapping in 
the sense that all maximally extended null-geodesics approach S- at one 
end and S+ at the other. The full set of geometric hypotheses is described 
in detail in Section 13.21 
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Figure 1. The polar and equatorial regions in Minkowski space 

We consider solutions w to the wave equation 

Dw = f€ C C °°(M°) 

on such a manifold so that w is tempered and vanishes near the "past in- 
finity" C- . In [20\ Section 5] the asymptotic behavior of the solution of the 
wave equation was analyzed in C+ on Minkowski space in a manner that 
extends to our more general setting in a straightforward manner, giving a 
poly homogeneous asymptotic expansion in the boundary defining function p; 
the exponents arising in this expansion are related to the resonances of the 
Laplace operator associated to a certain natural asymptotically hyperbolic 
metric on C+. 

The main result of this paper is to obtain the precise asymptotic behavior 
of the solution w at the light cone at infinity, S+ = <9C+, performing a 
uniform (indeed, conormal, on an appropriately resolved space) analysis 
as S+ is approached in different ways. This amounts to a blow-up of <S+ 
in M. In Minkowski space (t,x) € M 1+3 , locally near the interior of this 
front face (denoted ff), the blow up amounts to introducing new coordinates 
p = (r 2 + t 2 + l) -1 / 2 , s = t — r, y = x/r, the front face itself being given 
by p = 0, so s = t — r, y = x/r are the coordinates on the front face. 
More generally, if p is a defining function for the boundary at infinity of M 
and v is a defining function for S+, with (v,y) a coordinate system on dM, 
we can let s = v/p and use s,y as coordinates on the interior of the front 
face of the blow-up. Thus, s measures the angle of approach to S+, with 
s — > +oo corresponding to approach from C+, while s — > — oo corresponding 
to approach from Co- 

In order to make a statement without compound asymptotics, we consider 
the so-called radiation fields. Thus, in this paper we show the existence of 
the Friedlander radiation field 

TZ+[w](s,y) 

given by restricting an appropriate rescaling of the derivative of w to the 
new face obtained by blowing up the future light cone at infinity 5+; see 



Section 10 for the detailed definition. The function 1Z+ thus measures the 
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Figure 2. The radiation field blow-up of Minkowski space 

radiation pattern seen by an observer far from an interaction region; in the 
case of static metrics, it is known to be an explicit realization of the Lax- 
Phillips translation representation as well as a geometric generalization of 
the Radon transform [3J. 

Our main theorem concerns the asymptotics of the radiation field as s, 
the "time-delay" parameter, tends toward infinity, and more generally the 
multiple asymptotics of the solution near the forward light cone. 

Theorem 1.1. If(M,g) is a compact manifold with boundary with a non- 



trapping Lorentzian scattering metric as defined in Section 3.2 and w is a 
tempered solution of\D g w = f £ C^°(M°) vanishing in a neighborhood ofC-, 
then the radiation field of w has an asymptotic expansion of the following 
form as s —> oo: 

K+[w](s,y)~J2 E ^ K s-^-\\og S r 

j K<rrij 

Moreover, w has a full asymptotic expansion at all boundary face^\ with the 
compound asymptotics given by: 

N 

(1.1) !»~P^^^ Yl P%gp\ K (log p + log s) a a jeKa s- 1 ^ 

j £=0 K+a<rhjf 



Remark 1.2. Although it may appear in (1.1) that the log terms prevent 
the definition of the radiation field, we show in Section [8] that the log terms 
cancel in the £ = term, enabling the restriction to p = 0. 

Remark 1.3. In Minkowski space, the requirement that w vanish in a neigh- 
borhood of C_ implies that w is the forward solution of O g w = f. One 
should then think of the vanishing requirement as analogous to taking the 
forward solution of O g w = f. 

lr The power of s in this formula differs from the previous one by 1 due to a derivative 
in the definition of the radiation field. 
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A question of considerable interest is, of course, whether the radiation field 
actually decays as s — > +00 and, more generally, the description of the expo- 
nents <jj. Remarkably, these are the resonance pole^Jof a naturally-defined 
asymptotically hyperbolic metric defined on C+. We denote the family of 
asymptotically hyperbolic operators by ) and record the following corol- 
lary: 

Corollary 1.4. If there exists C > such that L~\ ^ as no P°l es a t °~ 
with Ima > —C then for all e > 0, the radiation field decays at a rate 
0{s- c ~ 1+ '). 

One class of spacetimes to which our theorem (and corollary) applies is 
that of normally short-range perturbations of Minkowski space, i.e., pertur- 
bations of the metric which are, relative to the original metric, 0(p 2 ) in 
the normal-to-the boundary component, 0(1) in the tangential-to-the- 

boundary components, ^5- and dv ^ v , and 0(p) in the mixed compo- 
nents. In particular, note that we are permitted make large perturbations 
of the spherical metric on the cap C+, hence in these "tangential" metric 
components our hypotheses allow a much wider range of geometries than 
even traditional "long-range" perturbations of Minkowski space. 

In the more restrictive setting of "normally very short range" perturba- 
tions, where we add an additional 0(p) of vanishing to all metric compo- 
nents, we recover the same asymptotically hyperbolic problem at infinity as 
in the Minkowski case, and therefore exhibit the same order decay as seen 
on Minkowski space. In particular, in odd spatial dimensions one has rapid 
decay of solutions of the wave equation away from the light cone. Thus, we 
obtain the following corollary for "normally very short range" perturbations 
of Minkowski space: 

Corollary 1.5. If (M,g) is a normally very short range non-trapping per- 
turbation oj 'n- dimensional Minkowski space, w vanishes nearC-, and\Z\ g w = 
f G C^°(M°), then the radiation field of w has an asymptotic expansion of 
the following form: 



K+[w](s,u>) 



0(s 00 ) n even 

,E^o£ K <i^t-J(logs) K a jK n odd 

More generally in the case of normally short-range perturbations, if the O(l) 
metric perturbations of the tangential-to-the-boundary metric components 
are sufficiently small then the radiation field still decays as s —> +00: 

^ + M( s ,^)< s - 2+e 



2 More precisely, the poles we are interested in are those of the inverses of a family of 
operators that looks to leading order like an asymptotically hyperbolic Laplacian. It is 
not in general a spectral family of the form P — a 2 however: the cr-dependence is more 
complex. 
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The polynomial decay of solutions of the wave equation may be com- 
pared with the Klainerman-Sobolev estimates [12] . (We refer the reader to 
the book of Alinhac [TJ Chapter 7] for a more detailed introduction to such 
estimates.) In re-dimensional spacetimes where the isometries (and confor- 
mal isometries) of Minkowski space (i.e., the translations, rotations, Lorentz 
boosts, and scaling) are "asymptotic isometries" (or "asymptotic conformal 
isometries" ) , then solutions w of the wave equation exhibit decay of the form 



In terms of these coordinates, the asymptotic expansion we obtain implies 
that on our class of Lorentzian manifolds (in particular, on perturbations of 
Minkowski space) , there is some a so that solutions w of the wave equation 
satisfy 



When there are no eigenvalues of the associated asymptotically hyperbolic 
problem, then a < 0, while on normally very short range perturbations of 
Minkowski space, a = —n/2 if n is odd and a = — oo if n is even. Further, the 
resonances of the asymptotically hyperbolic problem depend continuously on 
perturbations in an appropriate sense. The operator P~ l introduced below 
is stable, but may contain additional poles at certain pure imaginary integers 
as compared to the asymptotically hyperbolic problem (as is the case in 
even dimensional Minkowski space). Although such poles do not contribute 
to the asymptotics of the radiation field, under perturbations they may 
become poles of L~^. Thus, for small normally short range perturbations 
of Minkowski space, a is close to — min(2,n/2) (rather than — oo)j^| 

The class of spacetimes we consider is geometrically more general than the 
class of spacetimes on with the Klainerman-Sobolev estimates hold, but we 
require a complete asymptotic expansion of the metric (and thus consider- 
ably more smoothness at infinity). The methods we employ would, however, 
allow also for finite expansions when the metric has a finite expansion by 
more careful accounting. 

1.1. A sketch of the proof of Theorem |1.1| , We start with a tempered 
solution w of O g w = f G C£°(M ) vanishing identically in a neighborhood 
of C_ . We then fix \ £ C°° (M) supported near dM so that x is identically 
near C_ , identically 1 near the portion of the boundary where w is non- 
vanishing. In particular, the support of w d\ is compact in M° and x w = w 
near dM. We then consider the function u = ^f™ -2 )/ 2 ^ and set 



In higher dimensions, one may improve this statement to obtain a close to —n/2 by 
a careful analysis of resonant states supported exactly at the light cone. As the most 
interesting case is n — 4, when n/2 = 2, we do not pursue this improvement here. 





n b = P -V (n - 2)/2 n 9 p 



,(n-2)/2 
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The function u then solves D^u = f for some other function / £ C C °°(M°). A 
propagation of singularities argument (Section [4]) shows that u is conormal 
to {p = v = 0}. 

We now set P a = N(Oj,), where N is the reduced normal operator, i.e., the 
family of operators on the boundary at infinity obtained by Mellin transform 
in the normal variable. If we set u a , f a to be the Mellin transforms of u, 
and /, respectively, then u a solves 



PaU a = U+ errors, 



where the additional correction terms arise because is not assumed to be 
dilation- invariant. We show that the operator P a fits into the framework of 
Vasy [20j and modify the argument of that paper to show that P a is invertible 
on certain variable-order Sobolev spaces (Section [5]). The argument further 
shows that P~ l is a meromorphic family of Fredholm operators with finitely 
many poles in each horizontal strip. In fact, (most of) the poles of P" 1 
may be identified with resonances for an asymptotically hyperbolic problem 
(Section [7]). 

An argument of Haber-Vasy [6] implies that the residues of at the poles of 
P~ l are L 2 -based conormal distributions. In Sections [6] and [8] we show that 
they are in fact classical conormal distributions and thus have an expansion 
in terms of v. We calculate the leading terms of the expansion somewhat 
explicitly. Inverting the Mellin transform and iteratively shifting the contour 
of integration (Section [9| realizes these residues as the coefficients of an 
asymptotic expansion for u in terms of p. 

A slight complication is that not only the terms of the expansion become 
more singular as distributions on dM as one obtains more decay, which is 
indeed necessary for them to contribute to the radiation field in the same 
way (i.e. letting p — > with s = v/p fixed), but the remainder term also 
becomes more singular. We use the a priori conormal regularity, as shown in 
Section [4j to deal with this potential problem. The philosophy here is that 
since the algebra of b-pseudodifferential operators, discussed in Section [2] 
with further references given there, is not commutative to leading order in 
the sense of decay at dM (unlike, say, Melrose's scattering pseudodifferential 
algebra), one first should obtain regularity in the differential sense, which 
is the conormal regularity of Section |4j and then proceed to obtain decay, 
which is the ultimate result of this paper. 

Finally, rewriting the expansion in terms of the radiation field blow-up 
s = v/p yields an asymptotic expansion at all boundary hypersurfaces. The 
explicit computation of the leading terms shows that the logarithmic terms 
match up and thus u may be restricted to the front face of the blow-up, 
yielding the radiation field (after differentiation), and proving Theorem |1.1| 
in Section [TUl 
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2. B-GEOMETRY AND THE MELLIN TRANSFORM 

2.1. Introduction to b-geometry. For a more thorough discussion of b- 
pseudodifferential operators and b-geometry, we refer the reader to Chapter 
4 of Melrose [13] . 

In this section and the following, we initially take M to be a manifold with 
boundary with coordinates (p, y) G [0, 1) x K n_1 yielding a product decompo- 
sition M D U ~ [0, 1) X dM of a collar neighborhood of dM. In particular, 
for now we lump the v variable in with the other boundary variables as it 
will not play a distinguished role. 

The space of b-vector fields, denoted V&(M) is the vector space of vector 
fields on M tangent to dM. In local coordinates (p, y) near dM, they 
are generated over C°°(M) by the vector fields pd p and d y . One easily 
verifies that V&(M) forms a Lie algebra. The set of b-differential operators, 
Diff£(M), is the universal enveloping algebra of this Lie algebra: it is the 
filtered algebra consisting of operators of the form 

(2.1) A= Yl a, j>a (p,y)(pD p yD2em&Z(M) 

\a\+j<m 

(locally near dM) with the coefficients aj a 6 C°°(M). We further define a 
bi-filtered algebra by setting 

Diff^(M) = p- 1 DiffC 1 (M). 

The b-pseudodifferential operators Vl/i!(M) are the "quantization" of this 
Lie algebra, formally consisting of operators of the form 

b{p,y,pD p ,Dy) 

with b(p, y, £, rj) a Kohn-Nirenberg symbol; likewise we let 

*^'(M) = p- l ^(M) 

and obtain a bi-graded algebra]^] 

The space V&(M) is in fact the space of sections of a smooth vector bundle 
over M, the b-tangent bundle, denoted h TM. The sections of this bundle are 
of course locally spanned by the vector fields xd x ,d y . The dual bundle to 
h TM is denoted b T*M and has sections locally spanned over C°°(M) by 
the one-forms dx/x,dy. We also employ the fiber compactification h T*M of 
b T*M, in which we radially compactify each fiber. A set of local coordinates 
on each fiber near {v = p = 0} is given by 

V = -, £ = -, 7) = -. 

Ill 

The symbols of operators in ^ft(M) are thus Kohn-Nirenberg symbols 
defined on b T*M. The principal symbol map, denoted o\,, maps (the classical 



4 The convention we use for the sign of the weight exponent I is the opposite of that 
employed in some other treatments; we have chosen this convention as differential order 
and the weight order behave similarly: the space increases if either one of these is increased. 
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subalgebra of) ty™ ,l (M) to p 1 times homogeneous functions of order m on 
b T*M. In the particular case of the subalgebra Diff™''(M), if A is given by 
(2.1 ) we have 

a h (p- l A) = p- 1 Yl Ha(p,y)eri a 

\a\+j<m 

where £,77 are "canonical" fiber coordinates on b T*M defined by specifying 
that the canonical one-form be 

dp dy 
£— +v ■ — . 
P P 

In addition to the principal symbol, which specifies high-frequency asymp- 
totics of an operator, we will employ the "normal operator" which measures 
the boundary asymptotics. For a b-differential operator given by (2.1), this 
is simply the dilation-invariant operator given by freezing the coefficients of 
pD p and D y at p = 0, hence 

N(A) = a ^(°< y)(P D p) iD y G Dif C([0, oo) x dM). 

\a\+j<m 

The Mellin conjugate of this operator is known as the "reduced normal 
operator" and is simply the family in a of operators on dM given by 

\a\+j<m 

This construction can be extended to b-pseudodifferential operators, but we 
will only require it in the differential setting here. 

Here and throughout this paper we now fix a "b-density," which is to say 
a density which near the boundary is of the form 

dp 



f(p,y) 



A dyi A • • • A dy n - 
P 



with / > everywhere. Let L^(M) denote the space of square integrable 
functions with respect to the b-density. We H™(M) denote the Sobolev 
space of order m relative to L^(M) corresponding to the algebras Diff™(M) 
and ^™(M). In other words, for m > 0, fixing A G ^^(ilf) elliptic, one has 
w e H™(M) if w e Ll(M) and Aw £ Ll(M); this is independent of the 
choice of the elliptic A. For m negative, the space is defined by dualization. 
For m a positive integer, one can instead give a characterization in terms of 
Diff^(M). Let H™' l (M) = p l H^(M) denoted the corresponding weighted 
spaces. 

We recall also that associated to the calculus Vl/? (M) is associated a notion 
of Sobolev wavefront set: WF™' 1 ^) C h S*M is defined only for w G Hr°°' 1 
(since ^^(M) is not commutative to leading order in the decay order); the 
definition is then a £ WF™' (w) if there is Q G ^°(M) elliptic at a such 
that Qw G H™ ,l (M), or equivalently if there is Q' G ^™' l (M) such that 
Q'w G Ll(M). We refer to [TUl Section 18.3] for a discussion of WF b from 
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a more classical perspective, and [Ml Section 3] for a general description of 
the wave front set in the setting of various pseudodifferential algebras; |21|, 
Sections 2 and 3] provide another discussion, including on the b-wave front 
set relative to spaces other than L 2 (M). 

2.2. Scattering geometry. We now turn to the analogous concepts of 
"scattering geometry" which will be less used in this paper but which are a 
useful motivation. For a full discussion of scattering geometry, we refer the 
reader to Melrose |15| . 

In analogy to the space of b- vector fields, we define scattering vector fields 
as V sc = pVb] that is to say, the vector fields when applied to p must return a 
smooth function divisible by p 2 . They are locally spanned by p 2 d p and pd y . 
They form sections of a bundle SC TM; the dual bundle, SC T*M has sections 
locally spanned by dp/p 2 , dyj p. As motivation for our discussions of the form 
of the "scattering metrics" below, we remark that if we radially compactify 
Euclidean space, the constant vector fields push forward to be scattering 
vector fields on the compactification, hence sections of the tensor square of 
SC T*M are the natural place for asymptotically Euclidean or Minkowskian 
metrics to live. 

The scattering differential operators are those of the form (near dM) 
Y, a^(p,y)(p 2 D p y(pD y r E Diff£(M). 

\a\+j<m 

Again, this space of operators can be microlocalized by introducing scatter- 
ing pseudodifferential operators which are formally objects given by 

b{p,y,p 2 D p ,pD y ) 

with b(p, y, £, 77) a Kohn-Nirenberg symbol on the bundle SC T*M. There are 
of course associated scales of Sobolev spaces, which we will not have occasion 
to use in this paper, as well as wavefront sets which are described in detail 
in p2]. 

2.3. Mellin transform. We first recall the definition of the Mellin trans- 
form on R_|_. For a smooth compactly supported function u on IR + , u a = 
/ °° p~ ta ~ l u(p) dp. Because u is compactly supported, u a is an entire func- 
tion of a which decays rapidly along each line of constant Imo". Will also 
use the notation 

M.u = u 

for the Mellin transform. 

The Mellin transform on M + is equivalent to the Fourier transform by 
the substitution x = log p. In particular, the Plancherel theorem guarantees 
that it extends to an isomorphism of Hilbert spaces 

L 2 (R + ;dp/p) -> L 2 (R), 

and, more generally, to an isomorphism with a weighted space, 

p 5 L 2 (R + ;dp/p) -> L 2 {{lma = -6}). 
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Moreover, the Mellin transform intertwines pd p with multiplication by la: 



The inverse Mellin transform is given by integrating u a p la along a horizontal 
line {Im<j = C}, provided this integral exists. 

Near the boundary, we use the boundary defining function p to obtain a 
local product decomposition: M = [0, e) p x dM. This local product decom- 
position allows us to define the Mellin transform for functions supported 
near dM via cut-off functions that are identically 1 for p < e/2. In what 
follows, this definition suffices, as we may always cut off the functions in 
which we are interested away from the boundary. Note that this definition 
of the Mellin transform depends both on the boundary defining function p 
and on the cut-off functions chosen, but this dependence will not make a 
difference in the sequel. 

We additionally recall the space of L 2 -based conormal distributions 1^ 
on the boundary. Here we finally split the boundary coordinates locally 
into (v, y) G R x W 1 ^ 2 rather than using y to denote all of them. For the 
hypersurface Y = {v = 0} C dM, u G I^{N*Y) means that u G H s (dM) 
and A\ . . . A^u G H s for all k and for all Aj G ^ l {dM) with principal symbol 
vanishing on N*Y. 

We now record some additional mapping properties of the Mellin trans- 
form: 

Definition 2.1. Let C u denote the halfspace Ima > —u and let %{C U ) denote 
holomorphic functions on this space. For a Frechet space J 7 , let 



denote the space of g„ holomorphic in C u taking values in T such that each 
seminorm 



is uniformly bounded in v' > —v. 

Note the choice of signs: as v increases, the halfspace gets larger. 

We in fact allow elements of %(C„) to take values in a -dependent Sobolev 
spaces, or rather Sobolev spaces with u-dependent norms. In particular, we 
allow values in the standard semiclassical Sobolev spaces HV 1 on a compact 
manifold (without boundary), with semiclassical parameter h = |cr| _1 . Re- 
call that these are the standard Sobolev spaces and up to the equivalence 
of norms, for h in compact subsets of (0, oo), the norm is just the standard 
H m norm, but the norm is /i-dependent: for non-negative integers m, in 

coordinates yj, locally the norm ||<7||f/™ is equivalent to \JY1 \\hD yj g\\^ 2 • 

We will require some more detailed information about mapping properties 
of the Mellin transform acting on b-Sobolev spaces. 



{pd p u) a = iau a 



^(C y )n(a)- fe L 00 L 2 (M;^) 
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Lemma 2.2. Let u G H™' (M) be supported in a collar neighborhood of 
dM. Then 

Mu G H(Q) n ( ( r) max(0 - m) L oo L 2 (M;iJ m (X)). 
If u G H™' l (M) is furthermore conormal to p = v = i/ten 

TWu G W(Q) n (<7)-°°l°°l 2 (]R; /( m )(iv*y)). 

T/ie inverse Mellin transform maps 

H{C V ) n (a)- ^ 00 !/ 2 ^; J«(iV*y))) 

into 

p^°H^([0,oc) p ;I^(N*Y)) 
which in turn is contained in 

p a-0 v s-l/2-0 L oo 

Proof. For m a positive integer, the first result follows since lying in H™' 1 
implies that 

dfflu G p l Ll 

for all \a\ + \/3\ < m hence 

d^d^Mu G H(Q) n L°°L 2 (R; L 2 (X)); 

the result for general m > follows by interpolation. For m < 0, choose a 
positive integer rh such that m + rh > 0; then it can be written as a finite 
sum of terms of the product of at most rh b- vector fields applied to elements 
u' of H™ +m ' l (M). Now, the Mellin transform of such u' lies in H(Ci) Pi 
L°°L 2 (M;# m+m '(X)) by the first part; d y and d v act as vector fields on X 
and thus would lead to the conclusion that u is in "H(Q)nL°°L 2 (IR; H m (X)) 
if only they appeared; however, pd p Mellin transforms to ta, and thus we 
may obtain up to rh factors of a as well, leading to the desired weight when 
m < is an integer; interpolation gives the weight (without a loss) for all 
m < 0. 

The proof of the second and third parts is similar; here we use Sobolev 
embedding, and the fact that regularity under pd p , vd v and d y intertwines 
under Mellin transform with regularity under a, vd v , and d y . □ 

We remark further that Mellin transform maps H™'\M) into 
H(Q) n (ay^L^L^R; H°°(X)). 

This map is not onto, as there is no iterated regularity under pd v built into 
the latter space. 
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3. Geometric set-up 

3.1. Minkowski metric. As a preliminary to our discussion of Lorentzian 
scattering metrics, we record the asymptotic behavior of the Minkowski 
space on W 1 , endowed with the Lorentzian metric with the mostly minus sign 
convention (here we are following the notation of [20 J). We take coordinates 
t, x±, . . . , x n -i, and set 

t = p _1 cos 9, 
Xj = p~ l Uj sin 9, 

with ui £ § n ~ 2 . The Minkowski metric is then 

r 2 7 2 ( cos ^ ^9 ■ ^d9\ 2 sr^ ( „ dp dO 
dt — > dx = = sm9 — — > — a;,- sniff- + uja cosd h 

^ V P P ) ^\ 3 P 2 3 P 

dp 2 d9 2 n [dp d9 d9 dp\ 2 

= cos - cos 26—^- + sm20 ® 1 <8> - sur 

P A P 2 VP P P P J 

Here e?u; 2 represents the standard round metric on the sphere. 

As the function cos 29 clearly plays an important role here, we set 

v = cos 29, 

replacing the 9 coordinate by v, and write 

,„„. dp 2 v dv 2 1 ,dp dv dv dp s 1 — v dw 2 

3.1 g = V J^- + 

p 4 4(1 — f z j p z 2 v p z p p p A ' 2 p z 

We remark that this form of the metric in these extremely natural coordi- 
nates does not conform to the standard "scattering metric" hypotheses |15j 
often employed in the Riemannian signature, in which cross terms of the 
form (dp/ p 2 ) (8> (dy/p) with y a general smooth function are forbidden. 

3.2. General hypotheses. Let (M,g) be an n-dimensional manifold with 
boundary equipped with a Lorentzian metric g over M° such that g ex- 
tends to be a nondegenerate quadratic form on SC TM of signature (+, — 

Definition 3.1. We say that g is a Lorentzian scattering metric if g is a 
smooth Lorentzian signature symmetric bilinear form on SC TM, and there 
exist a boundary defining function p for M, and a function v G C°°(M) such 
that 

(1) with V a scattering normal vector field, g(V,V) has the same sign 
as v at p = 0, 

(2) in a neighborhood of = 0, p = 0} we have 

dp 2 ( dp a a dp. g 
p 4 p z p p p z ' p z 



^The hypotheses below imply that, even if the boundary is disconnected, v vanishes on 
each component of the boundary. 
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with a a smooth 1-form on M so that a|p=o is a 1-form on X = dM 
(i.e., the dp component of a vanishes at p = 0), and g a smooth 
symmetric 2-cotensor on M so that 

g\Ann{d P< dv) is positive definite. 
We further require that 

a = -dv + 0(v) + O(p) near v = p = 0. 

Remark 3.2. We remark that while it might be tempting to mandate also 
the vanishing of the dv 2 / p 2 component at v = as we have in the exact 
Minkowski case, this condition is highly non-invariant, in that it requires a 
product decomposition of X. 

Remark 3.3. The function v must have a non-degenerate 0-level set when 
restricted to X, since otherwise our metric would be degenerate at v = 0. 

Remark 3.4. Note that near v = 0, V = p 2 d p gives g(V,V) = v, which 
has the same sign as v, so the first and second parts of the definition are 
consistent. 

We make two additional global assumptions on the structure of our space- 
time: 

Definition 3.5. A Lorentzian scattering metric is non-trapping if 

(1) The set S = {v = 0, p = 0} C X splits into S+ and S-, each 
a disjoint union of connected components; we further assume that 
{v > 0} splits into components C± with S± = dC±. We denote by 
Co the subset of X where v < 0. 

(2) All null-geodesies tend to S± as their parameter tends to ±oo (or 
vice versa) j^] 

In particular, this implies the time-orientability of (M,g) by specifying 
the future light cone as the one from which the forward (in the sense of the 
Hamilton flow) bicharacteristics tend to S+. 

Remark 3.6. C stands for "cap" as in the Minkowski case C+ is simply the 
spherical cap \6\ < ir/4. The assumption the 5+ bounds a spherical cap 
is in fact not necessary for us to prove any of the Fredholm properties in 
however it is of course necessary to recognize the poles of the resulting 
operator as resonance poles on a cap, and hence in order to know that there 
are finitely many resonances in any horizontal strip in C, which is crucial to 
the development of our asymptotic expansion. 

One should think of the second assumption as a non-trapping assumption 
on the light rays in the spacetime. It is unclear whether the assumptions 



More precisely, we assume that all null-geodesics flow from the radial points at S- to 
those at S +, o r vice versa. Discussion of the radial points and their location is contained 
in Section 3.4 
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imply that forward tempered solutions, i.e. for compactly supported /, solu- 
tions of O g u = f with / = near C_, always exist, though it is not hard to 
show that there are only finite dimensional obstructions to solvability and 
uniqueness in fixed weighted spaces. 

Near v = 0, which is away from the critical points of v\x, we may choose 
yi, . . . ,y n -2 £ C°°(M) so that (v,y) constitute a coordinate system on X 
and (p, v, y) thus give coordinates on M in a neighborhood of X. Moreover, 
(p,v,y) also provide a product decomposition of that neighborhood into 
[0, e) p x X. In the frame 

p 2 d p , pd v , pd y , 

associated to these coordinates, the metric (when restricted to the boundary 
{p = 0}) thus has the block form 

/ v -^ + a v a\v ... a n ^ 2 v \ 
-\ + aov b c\ 
a±v ci — h\ t \ 



(3.2) 



Go 



ln-2V 

C n -2 

hn-2,1 



Cn-2 —hl : n-2 ■■■ —h n -2,n-2/ 

1) block negative definite, hence hij is positive 



\ a n -2V 

with the lower (n — 1) x (n 
definite. 

Blockwise inversion shows that in the frame 

dp dv dy 

P 2 P P 

the inverse metric when restricted to the boundary has the block fornj^] 

/ —q -2 + av -2T T + 0(v 

G 1 = -2 + av -4u + flv 2 -AvT T + 0(v 2 ^ 
\-2T + 0(v) -AvT + 0{v 2 ) -h- x + 0{y) 

In the above, h^ 1 = h l i is the inverse matrix of h™, q, a, (3, and T, are 
smooth near v = p = 0, and A T denotes the transpose of the matrix A. 

In a neighborhood of the boundary, i.e., at p ^ 0, there are further cor- 
rection terms in the inverse metric as the actual metric is given by 

G = Go + H, 

(0(p 2 ) O(p) 0(p)\ 
H = O(p) O(p) O(p) . 
\0(p) O(p) 0(p)J 

Thus in the inverse frame above, 



(3.3) 



G~ 



G o 1 



(0{p) O(p) 0(pf 

+ \0(p) 0(p 2 ) + 0(pv) O(p) 
\0(p) O(p) 0(p) / 
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The a here is a function and should not be confused with the 1-form a in the definition 
of the metric. 
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Abusing notation so that the above now denote functions of p as well as of 
v and y, in the actual coordinate frame d p , d v , d y , the dual metric becomes 

/gPPp 4 + 0{p 5 ) gP v p 3 + 0{p A ) gpyp 3 + 0(p 4 ) 

(3.4) gP v p 3 + 0(/j 4 ) g vv p 2 + 0(p 4 ) + 0(p 3 v) g^p 2 + 0(p 3 ) 

\g py p 3 + 0(p 4 ) g vy p 2 + 0(p 3 ) gyyp 2 + 0(p 3 ) 

where g" are given as above: 

g p P = -q g pv = -2 + av g pv = -2T + 0(v) 

' ' > " ) ' g m = -4v + (3v 2 g vy = -4vT + 0(v 2 ) g w = -h' 1 - 0(v) 

Again all terms are smooth. 



Cofactor expansion of equation (3.2) scaled to the frame d p , d v , d y shows 
that the determinant of the metric is 

\ g \ = p~ 2 ^)\G\ = p- 2 ^ {(f 2 - qv)\h\ + O(p)) 



2 d p \og\g\=-{n+l)p' 1 +0{l) 



In particular, 

1 

2 

l -d v \og\g\ = 0(1) 
l 2 d y \og\g\=0{l). 

3.2.1. Induced metrics. In this section we describe induced metrics on the 
"caps" C± (the components of {v > 0} bounded by S±) and on the "side" 
C ({« < 0}). 

We define the metric K on T* X via the inclusion r* : T*X ^ h T* x M 
(which is dual to the restriction map r : h TxM — > TX). As p 2 g is a b- 
metric, we define for to, rj € T*(X) the dual metric K^ 1 by 

K- 1 (u } ,r,) = -(p 2 g)- 1 (r*u } ,r* V )\ p= o. 

Observe that K~ l is the restriction of —(p 2 g)^ 1 to the annihilator of pd p 
(the b-normal vector field) at p = 0. 

The components of the dual metric K^ 1 are given in the frame d v , d y by 

j£W J£vy \ / _ vv _gvy 

K vy K yy ) = \ -g v y - 9yy 

where g" are the components of the dual metric of g in the frame p 2 d p , pd v , 
and pd y . 

Because p 2 d p is time-like near C± and K^ 1 is the restriction of —(p 2 g)^ 1 
to the annihilator of pd p , K^ 1 is nondegenerate, and, in fact, Riemannian 
in C±. In coordinates (v,y), the metric K on TX is given by 



I n— 2 n— 1 

K = — (1 + O(u)) cfo 2 +J^ 0(1) (dw (8) c% + dy + j (8) cfc)+ ^ K^dy^dyj. 



3=1 ».J=1 
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It is easy to see from the above expression that the metric 

k± = -K\c± 
v 

is an asymptotically hyperbolic metric (in the sense of Vasy [20]) on C±. 
Setting v = x 2 in this region ensures that k is an asymptotically hyperbolic 
metric (in the sense of Mazzeo-Melrose [13]) which is even in its boundary 
defining function (cf. the work of Guillarmou [5]). 

Similarly, because p 2 d p is space-like near Co, K~ l \c is Lorentzian (with 
the "mostly-plus" convention), and 

k = -K\ Co 

v 

is an even asymptotically de Sitter metric (with the "mostly-minus" conven- 
tion, as v < here) on Co- Indeed, if v = —x 2 , then the metric has the form 
used by Vasy [22j. The non-trapping assumption (2) above implies that the 
metric satisfies the conditions in Vasy's definition of an asymptotically de 
Sitter metric. 

The p components of the dual metric of g are also related to the com- 
ponents of the dual metric of K. In the p 2 d p , pd v , pd y frame for g and the 
d v ,dy frame for K, we have 



g PP = - l -(AqK™ + 0{v 2 )), g Pv = -±-(K™ + 0(v 2 )), g"> = -i- [IT* + 0(v 2 )) 

As K~ l is the lower-right block of —g~ l and g pp = v, the volume forms 
of g and K (and hence the asymptotically hyperbolic and asymptotically de 
Sitter metrics k±, ko are also related: 



= ( v^y/m + O(P)) = p-^ V n ' 2 J\k ± ,o\ + 0{p) 



3.3. The form of the d'Alembertian. In this section we compute the 
form of the operator D g and its normal operator N(p~ 2 D g ). 

Putting the calculations of the metric components and the volume form 



in Section 3.2 together, we compute the form of □ near p = (here we use 
VG = p n+1 y/g and recall that g" be given by (3.5)): 



{gPP + O(p)) (pd p ) 2 + + 0{p)) (pd p )d v + (g^ + 0{p)) {pd p )d y 
+ (2 - n) ((gPf + O(p)) pd p + QT + O(p)) d v + (g^ + O(p)) d y ) 
+ -±=d v [(g pv + 0(p)) VGpd p + (g vv + O(p)) VGd v + (g^ + O(p)) VGd, 



+ -La y (y» + o( P )) VGpd p + (g v y + o( P )) Vgo v + ( g *y + o( P )) VGd y ) 
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Adopting now the notation of Vasy [20 , 
-P. = -N(p- 2 O a ) 

1 <8JrVG8 v+a -»VG8 v ) + B J g ^ + ^V08, 



Vg l 

+ g pv (2ia + 2-n)d v + g py (2ia + 2-n)d y 
+ 10 -±=d v (g pv Vo) + ^=dy (g py Vc) + g pp ia 
In particular, near v = 0, 

-p a = (-4« + o{v 2 )) d 2 + o{v)d v d y - (h ij + o(v)) d yr dy 3 + o(i)d y 

+ 2 (n - 4 - 2ia + O(u)) d„ + g(<r), 

with g a smooth function in u and y with values in quadratic polynomials 
in a. 

In our asymptotic expansions (and in the analysis of the radiation field), 
it is more convenient to deal with 

(3.6) p a = n ( P -("- 2 )/V 2 n s p ( "- 2)/2> 

than with P a , in part to more directly correspond to the setting of [20]. To 
this end we note simply that 

(3-7) P a = PcT-i(n-2)/2 

hence since the d 2 and d v terms of P c may be written 

(3.8) - 4 (jv + 0(v 2 ))D 2 v + Q(n - 4 - 2ia) + 0{v)) 

we have 
(3.9) 

P(j = -4 ((v + 0(v 2 )) D 2 V + ({a - i) + 0(v)) D v )+0(l)d 2 y +0(l)d y +0{v)d v d y +0{a 2 ). 

3.3.1. Relationship with the induced metrics. In the regions C± and Co of 
the boundary, P a may be written in terms of the metrics k± and ko. 

We first work near C±. By an explicit computation, there is a (<r- 
dependent) vector field X(a) tangent to v = and a (cr-dependent) smooth 
potential V(a) E C°°{X) so thafl 

1 n I itr n itj 1 { n (ft — 2)^\ , , 

(3.10) v ^v* + ^P ff v~*~^v' = -A fc± + ( a 2 + - — J +vX(a) +vV{a). 

(In terms of the variable x given by v = x 2 , the vector field X is in fact a 
0- vector field in the sense of Mazzeo-Melrose |13j.) Moreover, if all a,j and 

The result of this computation should perhaps not be too surprising, as the entries 
of the inverse metric of k agree up to a factor of v with a block of the inverse metric of 
g, accounting for the second-order terms. Moreover, it is easy to check that the operator 
on the left side is a b-differential operator. The remainder of the computation requires 
checking only that the b-normal operators of the two sides agree. A similar computation 
is carried out in |201 Section 5]. 
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q vanish identically on X (as is the case in Minkowski space and normally 
very short range perturbations of Minkowski space) then X = and V = 0. 
We now consider Cq. The same calculation as above implies that 

(3.11) \ v \*+*+-* P„\v\*~*~*- = Dfeo + ( ^ + 4 J +u^(a)+uF(c7), 

where X and 1/ are as above. In particular, P a is a hyperbolic operator on 
Co and an elliptic operator on C±. 

3.4. Location of radial points. We now study the flow associated to the 
Hamilton vector field of P a . In particular, we are interested in the radial 
points of the vector field, i.e., those points in the characteristic set where it 
is proportional to the fiber-radial vector field. As P a is hyperbolic for v < 
and elliptic for v > 0, the only possible radial points must occur when v > 0. 

As is not a critical point of v, we may take 

^ydv + r] ■ dy 

to be the canonical one- form on T*X, the principal symbol of P a is 

a(P a ) = -(4v - f3v 2 ) 7 2 - (4vT + 0(v 2 ))jr] - (h» + Olv))^. 
Letting H denote the resulting Hamilton vector field on T*X, we have 

(3.12) = (-4vy + /3 2 7 + vr] ■ T)d v + (vyT, + g y ^)d y] + .d 7 + 

with the • terms homogeneous of degree 2 in the fiber variables. We now 
analyze the radial points of the vector field. The components in the base 
variables are given by 

(-4vy + /3v 2 7 + vr) ■ T)d v + (vjTj + g Vi ^)d yj . 

These coefficients must vanish at the radial set, which we have already ob- 
served to lie over v = 0. In particular, we must have 

g ViV iTh = 

for all j. As g Viyj is nondegenerate at v = 0, we must have rj = on the 
radial set. 

We now easily verify that indeed the vector field at points 

V = 0, T] = 

is radial; hence these are in fact the radial points. 

3.5. Structure near radial points. We now verify several of the hypothe- 
ses of [20] near the radial points under the foregoing hypotheses. We have 
established that the radial points occur at 

L± = {v = 0, tj = 0, ±7 > 0} C T*X. 

We must now verify the following: 
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(1) For a degree —1 defining function p^ of S*X inside the fiber-radial- 
compactification of T*X, we have 

PooU Poo \ L± = =FA), Po € C°°(L±), fa > 0. 

(Equation (2.3) of [20].) 

(2) There exists a non-negative homogeneous degree function p$ van- 
ishing quadratically and non-degenerately exactly at L± and a /3i > 
such that 

^F/QooH/Oo — /3iPo ^ modulo cubic terms vanishing at L± 
(Equation (2.4) of [20].) 



To deal with the first property, we remark that from (3.12), we have 
= (2 7 2 +0(i>)+0(r ? ))a 7 +(0(7 ? 2 )+OM+0(i> 2 ))^+(-4z;7+/37; 2 7+^-T)a t ,+.a y 

where the big-Oh terms all have the homogeneities in 7, r\ required to make 
the overall vector field homogeneous of degree 1. Near 77 = we may employ 
the homogeneous coordinates 

1 AT V 

Poo = ri> N = n 

on the radial compactification of T*X, hence we compute that near L± 

l -H= p-J{{^2 + Oiv) + 0{N)) Poo d Poc 
( 3 - 13 ) + (t2N + 0(v 2 ) + 0(vN) + 0(N 2 ))d N 

+ (t4±Pv + n -r)vd v + »d y ), 



\L± 



hence 

=FPooHp c 
i.e., the first property holds with 

A) = 4. 

To verify the second property, we take, in our compactified coordinates, 

po = v 2 + N 2 . 



Applying (3.13) yields 



PooH(po) = =F(16u 2 + 8iV 2 ) + cubic terms in (v, N), 



i.e., 



= FPooH(/9o) — 8po = cubic terms in (v, N), 
hence the second property is satisfied with fi\ = 8. 
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We may thus compute the subprincipal symbol of P a (and hence of P a 
in terms of /3q. Indeed, we compute 

(3.14) = (2i)- 1 ((8 + i(n-4 + 2Ima) + 0(v))d v + 0(l)d y ) + O(l) 

= -2i((n - 2 + 2Ima) + 0{v))d v + 0(l)d y + 0(1) 
and consequently 

*{fy)-\P a - P;))\ v=0>v=0 = ±4 f-^- 2 - Ima) | 7 | 

(3 ' 15) / n-2 \ 

= ±M z--Tm* | 7 |. 



Note that, even apart from the shift by (n — 2)/2, the sign of (2i)~ 1 (P CT — P*) 
is switched as compared to [20] (where ±/3o Im cr|^y| was used with the present 
notation here). Switching the roles of P a and P* reverses this sign, and 
thus what we do here corresponds to what was discussed in [20] for the 
adjoint operator in the context of the general theory, though this reversal 
was pointed out there already in the Minkowski context in Section 5. 

Returning to the operator P a itself, we compute for later reference that 
by ( [335] ) and gjft , 

(3.16) ^ = ±^0^-^)1^ 

( n ~ 2 ) _ Im ( a _ t ( n -2)/2))=- Im a. 



3.6. b-radial points. It is also useful to compute the full b-structure of 
the radial set of 

L = p -(n-2)/2 p -2 Dp (n-2)/2 £ Diff 2 (M) _ 

in h T*M . Note that the powers are chosen here so that L is formally self- 
adjoint with respect to the b-density 

P n \dg\. 

The b-principal symbol of L is the 'same' as the sc-principal symbol of □ 
under the identification of h T*M and SC T*M, namely it can be computed 
to give 
(3.17) 

A = a b (L) = go?? - (4« - fiv 2 + 0{pv) + 0(p 2 )) 7 2 - 2(2 - av + 0(j>))fr 
+ 2gpy ■ rit + (2vT + O(p)) ■ 177 + //" " //,//,. 

where we write b-covectors as 

j. dp 

4 h 7 av + 77 ay. 
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The b-Hamilton vector field of a symbol A is 

(3.18) (<%A)( P cg + (d 7 \)d v + (d v x)d y - ( P d p x)d ( - (cU)<9 7 - (^a)^, 

so in our case we obtain 

H A =(2gP»£ + 2g™r, - 2 7 (2 - av + 0(p)))(pd p ) 

- 2((4« - (3v 2 + 0{pv) + 0(p 2 )) 7 

+ (2 - av + 0(p))£ + (vT + 0{p))rj)d v 
+ 2{gPy^ + (vT + 0(p)) 7 + g yiVj Vj)d y 

- (pd p X)d^ - {d v X)dj - (dyX)d v . 

Notice that this is a vector field tangent to Tq M M, with a vanishing <9,t 
component at p = 0; thus at h Tg M M this is precisely the same result as 



in (3.12) if one includes £ as a large parameter there, i.e. if one performs 



the semiclassical calculation. (The stated version of (3.12) corresponds to 
letting £ = 0.) In particular, the radial set 1Z within p = is exactly 
v = 0, r] = 0, £ = 0. Further, there are no radial points in p > 0, since 
the metric is a standard Lorentzian metric there (and there is no distinction 
between b-metrics and standard metrics in the interior). Now, on the fiber 
compactification of h T*M near 1Z we can use local coordinates, 

1 c Z - V 

v = -, £ = -, r/ = 

7 7 7 

to obtain the linearization of H A . That is, £ Vb( b T*M), i.e. is tangent 
to both p = 0, defining <9M, and = 0, defining fiber infinity, vanishes 
at dlZ (fiber infinity of the radial set), thus maps the ideal X of C°° func- 
tions vanishing at a point q G <97£ to themselves, and thus X 2 to X 2 , so it 
acts on X/X 2 = T* h T*M. In computing this, terms of vV\\ which vanish 
quadratically at dlZ can be neglected; modulo these we have 

uH x = -Apd p + (-8v - Ai)d v + 2(gpyi + v T + pc + g ViV >f)j)d yi 

- A(ud v + + fjdfj) + 1 2 V(^M), 

with c smooth. Correspondingly, the eigenvectors and eigenvalues of z^H^ 
are 

dv + d£, with eigenvalue — 8, 

(3.19) dp,du,d^,dr], with eigenvalue —4, 

2 dy + g Pi/ d£ + T du — c dp + dr/j , with eigenvalue 0. 

3.7. The radiation field blow-up. Although the existence of the radia- 
tion field for tempered solutions of D g w = f £ C^°(M°) with appropriate 
support properties is a consequence of our theorem, in this section we recall 



the definition of the radiation field for metrics of the form in Section 3.S ' 



^In less general settings, one may take the appropriate term of the expansion we obtain 
in the theorem as the definition of the radiation field. 
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We also assume a support condition on the solution w analogous to the 
one satisfied by the forward fundamental solution in more familiar contexts. 
More specifically, we suppose that 

(1) The metric g is a symmetric section of ( SC T* M)® 2 with signature 
(+,-...,-), 

(2) is not a critical value of v = g(p 2 d p , p 2 d p ), and in a neighborhood 
of v = 0, g can be written in the form 

dp 2 f dp a a dp\ q 

P 4 VP P P P J P 

(3) a = Tjdv + 0(v) + 0{p) near 5, i.e. near v = p = oj^j 

(4) (M°,g) has non-trapped null-geodesics (in the sense of the previous 
section). 

(5) The function w solves D g w = f 6 C^°(M°) and there is an sq so 
that near S+, w vanishes identically for v/p > sq. 

Assumptions (1)— (4) imply that S = {v = p = 0} is a smooth submanifold 
of M, while the support assumption (5) is analogous to requiring that w be 
the forward solution of O g w = f. 

We now blow up S = {v = p = 0} by replacing it with its inward 
pointing spherical normal bundle. This process replaces M with a new 
manifold M = [M; S] on which polar coordinates around the submanifold 
are smooth, and depends only on S (not the actual functions v and p). The 
blow-up comes equipped with a natural blow-down map M — > M which is 
a diffeomorphism on the interior. M is a manifold with corners with two 
boundary hypersurfaces: bf, the closure of the lift of X \S to M; and ff, the 
lift of S to M. Further, the fibers of ff over the base, S, are diffeomorphic 
to intervals, and indeed, the interior of the fibers is naturally an affine space 
(i.e. these interiors have M acting by translations, but there is no natural 
origin) . 

Given v and p, and not just S (but local coordinates on S are still irrele- 
vant), the fibers of the interior of ff in [M; S] can be identified with M, via 
the coordinate s = v/p. In particular, d s is a well-defined vector field on the 
fibers. 

We define "polar coordinates" 

R=(v 2 + p 2 ) 1/2 G[0,oo), 6 = ^^!, 

which are smooth on M. Near the interior of ff, we use the projective 
coordinates p, s = - p as well as local coordinates y on S. In these coordinates, 
a simple computation shows that the unbounded terms of p 2 g cancel near 
p = and hence p 2 g is a smooth Lorentzian metric in a neighborhood of the 
interior of ff (i.e., down to p = 0). 



If a is a different (but still constant) multiple of dv at p — v — 0, in some cases 
one can still perform an inhomogeneous blow-up (with inhomogeneity determined by the 
constant) to define the radiation field. 
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Given a solution w(p,v,y) of O g w = f with / smooth and compactly 
supported, we define the function 

u(p, s, y) = p~^w(p, ps, y). 

The wave operators for the metrics g and p 2 g are related by the somewhat 
remarkable identity 

2 — n 2— n / n — 2 \ o { "+ 2 2 — n \ 

P 2 DjW = P 2 Dp 2 Uj = p U p 2 g U - [p 2 U p 2 gP 2 j U 

we refer the reader to [1] for the details of this computation. Note that 7 is 
smooth on M because p is. Moreover, p 2 g is a nondegenerate metric near 
the interior of ff , and so O p 2 g — 7 is a nondegenerate hyperbolic operator 
near ff . This calculation thus shows that if w is a solution of Uw = f smooth 
compactly supported /, vanishing identically for s < sq, then the argument 
of Friedlander [4, Section 1] shows that w may be smoothly extended across 
ff . In particular, w and its derivatives may be restricted to ff . Note that the 
condition on the support of w is analogous to the support condition satisfied 
by forward solutions of the inhomogeneous equation. (The argument applies 
equally well to solutions of the homogeneous initial value problem with the 
same support property on globally hyperbolic spacetimes of this form.) 

Definition 3.7. If w is a solution of O g w = f satisfying the above support 
property, with / smooth and compactly supported, we define the (forward) 
radiation field of w by 

■ft+M( s >y) = d s u(0,s,y). 

Remark 3.8. Note that the smooth expansion of w across ff implies that it 
does not have singularities at s = 0. 

4. Propagation of b-regularity 

In this section we prove an initial conormal estimate for tempered solu- 
tions w of Ow = f G C C °°(M°) vanishing near C_. This estimate is used to 
begin the iterative scheme in Section |9j 

The basic background in this section is the propagation of b-regularity 
away from radial points (see, e.g., [H]), which we briefly recall here. Let 
L G ^ r (M), and let E C b S*M denote the characteristic set of L. 

Proposition 4.1. Suppose w G Hr°°'\M). Then 

(1) Elliptic regularity holds away from E, i.e., 

WF™V) C WF™- s > l ~ r (Lw) U S, 

(2) In E, WF™' ! («j) \ WF™~ s+1,r ~ l (Lw) is a union of maximally ex- 
tended bicharacteristics, i.e., integral curves of 
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Note that the order in WF™~ s+1,r ~ l (Lw) is shifted by 1 relative to the 
elliptic estimates, corresponding to the usual hyperbolic loss. This arises 
naturally in the positive commutator estimates used to prove such hyperbolic 
estimates: commutators in ^f^(M) are one order lower than products in 
the differentiability sense, hence the change in the first order relative to 
ellipticity, but not in the decay order, so there is no change there (cf. the 
radial point estimates below). 



We now turn to the radial set TZ, where Proposition 4.1 does not yield 
any interesting statements, and more refined arguments are needed. 

Definition 4.2. Let M C ^\{M) denote the *°(M)-module of pseudodiffer- 
ential operators with principal symbol vanishing on the radial set TZ = {p = 
0,^ = 0,^ = 0,^ = 0}. 

Note that a set of generators for M. over ^g(M) is pd p , pd v , vd v , d y (with 
symbols £, /r/, vy, m 7 enters to convert p and v to first order operators) and 
Id. 

Lemma 4.3. The module Ai is closed under commutators. 



Proof. While this can be checked directly from (3.18), a more conceptual 
proof is as follows. 

Recall that in the setting of a manifold without boundary, Ai being de- 
fined by the principal symbol of its elements vanishing on a coisotropic (i.e. 
involutive) submanifold would be closed under commutators. Here we are in 
an analogous situation, but need to employ the language of contact geome- 
try instead. We remark that while the symplectic form on h T*M is singular 
at the boundary, SC T*M is equipped with a smooth symplectic form sc wJ^] 
The submanifold h Tg M M is a contact manifold with contact form 

a = d(, — 7] ■ dy — 7 dv, 

given by contracting sc u; with p 2 d p and then multiplying the result by p. 
Although this is not a natural object on b T*, it is fixed once p is fixed. 
The principal symbols of elements of Ai are exactly those functions vanish- 
ing on the Legendrian submanifold TZ. We further note that the Hamilton 
vector field (w.r.t. the symplectic structure on the interior of h T*M) of a 
function vanishing on such a Legendrian submanifold is then tangent to it. 
Commutators of such vector fields then preserve the tangency. □ 
Proposition 4.4. Let L = p -(™-2)/2 p -2 Dp (n-2)/2 £ Diffg(M). If w G 
H™ ' (M) for some mo, I with mo + I < 1/2, Lw G H^~ 1,1 and w G H™' 1 
on a punctured neighborhood U \ dTZ of dlZ in h S*M (i.e. WF™' (kj) D (Z7\ 
&R) = %) then for m' < m with m' + I < 1/2, w G H™'' l (M) at dTZ (i.e. 
WF™''(to) n dTZ = 0,), and forN^N with m! + N < m and for A G M N , 
Aw is in H™'' l (M) at dTZ (i.e. WF™'' l (Aw) C\dTZ = $). 



we have fixed a defining function p on which the metric g depends, the "scattering" 
objects are in some sense more natural in our context. 
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Remark 4.5. In the situation that we care about, 7Z = 7Z+ U 1Z- splits into 
two components ("future" and "past") and we note that the proof in fact 
shows that the result holds at each component separately. 

This result is analogous to [B] , except p = produces an extra boundary 
(so we are in codimension 2), and 1Z is not Lagrangian ( b T*M is not sym- 
plectic at the boundary). The relevant input of the Lagrangian nature in [BJ 
is the eigenvectors and eigenvalues of the linearization, hence much the same 
proof goes through. This is also analogous to the 'easy' part, Section 11, of 
[IB] , describing the propagation of edge singularities, except here we have a 
source/sink rather than a saddle point, and thus the treatment is simpler. 

Proof. First we ignore the module. One proceeds inductively to prove 
WF™' (w)ddTZ = assuming that one already has shown WF™ ' (w)ndlZ = 
with m" = rh — 1/2. One can start with rh = min(mo + 1/2, m'), and in- 
creasing m by < 1/2, one reaches rh = m! in finitely many steps. 
Thus, one considers A £ ^ r (M) = p~ r ^l(M). Then 

%[L,A] € *£ +1,r (M), a h {i[L,A]) = H A a, a = a h , r , s (A). 

We choose 

a = p~ r v~ s <j) 2 , 

where (ft > 0, <p = 1 near 7Z, supported in U (supp eft will be further con- 



strained below). By (3.19), vH\a = (4(r + s) + c)a + e, where c vanishes at 
v = 0, and e is supported in supp dep. We take r + s < 0, and we choose the 
support of 4> so that |c| < |r + s\ on the support of <j>. Note that r + s < 
means vW\a necessarily has negative sign at least in some place on supp d<fi, 
since <fi has to increase along the flow as it approaches 1Z. Then we have 
v\-\\a = -b 2 + e, with b elliptic near K. Then with B € ^ { h s+1)/2,r/2 (M) 
with principal symbol b and WF' h (B) C supp b n h S*M (so for instance B 
can be a quantization of 6), and similarly with E £ ^ f ^ +1 ' r (M), 

i[L, A] = —B*B + E + F, F £ 

This gives an estimate 

\\Bwf < \(Ew,w)\ + \(Fw,w}\ + 2\(Lw,Aw}\ 

when w is a priori sufficiently regular. Given rh, I, we now take s = 2rh — 1, 
r = 21, so s + r < indeed. Note that F has order < 2m" , so the inductive 
assumption gives a bound for | {Fw, w) |. A standard regularization argument 
can be used to complete the proof: for instance one can use a regularizer 
ij)eiy) = (l + e^ -1 ) -1 = ^r^, e > 0, which is in S~ l for e > and is uniformly 
bounded in S° for e G (0, 1]. Thus, one lets 

then ^H;^ = ev~ 2 ^ 2 (vW\v") shows that the contribution of the regularizer 
to the principal symbol of the commutator is the negative of a square, pro- 
vided again that has sufficiently small support, i.e. adds another 'good 
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term' beside — b 2 . One can drop the corresponding term in the inequality 
given by quantized version, 

\\B e wf < \{E e w,w)\ + \(F e w,w)\ + 2\(Lw,A £ w}\, 

where the calculation (involving the pairing) now makes sense for e > 0. Now 
letting e — > the right hand side remains bounded, while B e —> B strongly 
in £(L^(M)), so one concludes Bw 6 L 2 (M) as usual. This completes the 
proof in the case when N = 0, i.e., when we have included no factors from 
the module M. in the test operator. 

In the general case one employs the methods developed by Hassell, Mel- 
rose and Vasy [8], adapted to a similar, but different (edge), setting by 
Melrose, Vasy and Wunsch in the appendix of [16]. For this purpose one 
uses generators of the module, denoted by Go = Id, Gi, . . . , G n , G n +i = AL, 
where A £ ^j^" 1 is elliptic near 1Z. A sufficient condition for these methods 
is that for i = 1, . . . , n, 

(4.1) iA[G i ,L]=J2C ij G j , 

j 

where 

(4.2) a bAQ (Cy) \n = 0. 

In our case this sufficient condition is satisfied by choosing dgi to be an 



eigenvector of v\\\ at 1Z, with eigenvalue —4 (cf. (3.19)), where Gi has 
principal symbol v~ 1 gi- Since dv and dgi have equal eigenvalues then, the 
conclusion for Cij follows. 

We thus prove iterative regularity under M. inductively in the power of 
the module as follows: we repeat the previous commutator argument, but 
with the commutant A replaced by 

Op(v / a)*(G a )*(G Q )Op( v / a) 

where G a = G" 1 . . . G^T^ 1 denotes a product of powers of the generators 
of Ai, hence G a G A^' a '. Considering all of these commutators at once, as 
G a runs over a basis of Ai N /.M^ -1 , we then follow the same argument 
as used when N = but now with systems of operators, taking values in 
with d = dim M N /M N ~ l . The main term in the commutator, arising 
from the commutators [L,Op(T/a)], is diagonal and positive, just as before 
(again, because the factor 4(r + s) + c is negative). Moreover the condition 



(4.2) permits us to absorb into this positive term those new terms that arise 



from commutators of L with G a and that have the maximum number of 

1 1 2 

module factors. Thus we are in the end able to estimate the terms ||-BG Q i(;|| 
(with B as before) where \a\ = N by terms microsupported away from 1Z 
and by terms involving G^w with \j3\ < N — 1, thus proving the result 
inductively. □ 



An immediate corollary is the following: 
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Corollary 4.6. Let L = p -{n-2)/2 p -2 Up {n-2)/2 £ Diffg(M), and let ir : 
b T*M — )■ M be the projection. Suppose w G HC°°' l (M) for some I, Lw G 
H™- L >\M). Suppose U is a neighborhood of ir{dlZ) and that all bicharac- 
teristics (in T,) of L that enter U, other than those in 1Z, possess a point 
disjoint from WF™''(w). Then for m! < m with m! + l < 1/2, w is in H™ ' 

on U and for A r G N with m' + N < m and for A G M N , Aw is in H™ ' l on 
U. 

Note that the hypotheses of the corollary at the future radial set hold 
automatically if L is non-trapping, i.e. all bicharacteristics tend to the future 
and past radial sets in the two directions of flow, and if w vanishes near S-. 

Remark 4.7. Corollary |4.6| implies that w is in fact conormal to the front 



face of the blow-up defined in Section 3.7 In particular, this implies that 
the coefficients in the asymptotic expansion of Theorem 1.1 may be taken 
to be smooth. 



5. The mapping properties of P a 

Having verified that the operator P a satisfies many of the hypotheses of 
the theorem of Vasy [20] , we now show that P a is Fredholm on appropriate 
function spaces. In this section we modify the argument of [20j to our current 
setting. 

With our global assumptions, the characteristic set of P a in S*X, X = 
dM, has two parts £± (each of which is a union of connected components) 
such that the integral curves of the Hamilton flow in £± tend to S± as the 
parameter tends to +oo. Writing the radial sets (for the Mellin transformed 
problem) at future, resp. past, infinity as A + , resp. A - (and within each one 
has two (unions of) components, = n S±), one is interested in the 
following two kinds of Fredholm problems, in which one requires a relatively 
high degree of regularity at A + , resp. A - , but allows very low regularity at 
the other radial set, A~, resp. A + . m 

To make this into a Fredholm problem it is convenient to introduce vari- 
able order Sobolev spaces and variable order ps.d.o's. This was originally 
done by Visik, Eskin [23], Unterberger [T7] and Duistermaat [2], and we 
recall this theory in Appendix lAl The main result that we use is Proposi- 



tion A.l, which shows that standard propagation of singularities arguments 



along forward null-bicharacteristics hold with respect to the spaces H s with 



We recall that in [20] such issues were avoided by using complex absorption arranged 
so that the resulting operator is elliptic at one of the radial sets, say A - , but is unchanged 
near A + . Thus, each bicharacteristic enters the complex absorption region in either the 
forward or backward direction, where the operator becomes elliptic due to the imaginary 
part of its principal symbol, hence only A + acts as a radial set for the operator with 
complex absorption added, and one could use standard Sobolev spaces as one did not 
have to deal with different regularity thresholds at A + and A~. 
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s £ C°°(S*X) defining the variable order, provided s is non-decreasing along 
the Hamilton flow. 

Now we recall, as computed in (3.16), that the quantity 

^( a )=±^a 1 (P a -P*)\ A ± 

is given the "constant" value — Ima : it is independent of the point in A 1 * 1 . 
Here the ± at the front of the right hand side corresponds to T,±, i.e. the 
subscript o/A^. Let 

-s ± {a) = \-P ± {a)= l -+\raa 

denote the threshold Sobolev exponents at A 1 * 1 , i.e. at the future and past 
radial sets. Thus, 

s+(cr) = s~(a), 

but this is actually not important below. Let Sft r be a function on S* X, 
such that 

(1) Sf tr is constant near A 1 * 1 , 

(2) Sf tr is decreasing along the H p -flow on £ + , increasing on XL, 

(3) Sft r is less than the threshold exponents at A + , towards which we 
propagate our estimates, i.e. Sf tr |A+ < s + (a), 

(4) Sf tr is greater than the threshold value at A - , away from which we 
propagate our estimates, i.e. Sf tr |A- > s~(a), 

then one can propagate regularity and estimates from A~ to A + as in [20} 
Section 2.4], and for all N (in practice taken very large) obtain estimates 

(5.1) \\U\\Htb < C(\\P a U\\ 

(More generally, the Sobolev exponent on the first term on right hand side 
would be Sft r — m + 1 where m is the order of P a \ here of course m = 2.) 
On the other hand, if s pas t is a function on S*X, such that 

(1) s past is constant near A 1 * 1 , 

(2) Spast is increasing along the H p -flow on E + , decreasing on E_, 

(3) Sp as t is less than the threshold exponents at A - , towards which we 
propagate our estimates, i.e. s pa st|A- < s~(cr), 

(4) Spast is greater than the threshold value at A + , away from which we 
propagate our estimates, i.e. s pas t|A+ > ^ + ( c7 )) 

then one can propagate regularity and estimates from A + to A - , and for all 
obtain estimates 

[|Z7||#«past < C (\\P a U\\ H 3 pllBt -, n +l + ||C/||^-jv). 

With s^'*{a) denoting the threshold Sobolev exponents for P*, the same 
considerations apply to P*, i.e., if s* ast is a function on S* X such that 

(1) s* ast is constant near A±, 

(2) s* ast is increasing along the H p -flow on S + , decreasing on S_, 
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(3) s* t is less than the threshold exponents at A , towards which we 
propagate our estimates, i.e. s* ast |^- < s~'*(a), 

(4) s past i s greater than the threshold value at A + , away from which we 
propagate our estimates, i.e. Sp ast |A+ > s + '*(cx), 

then one can propagate regularity and estimates from A + to A - , and for all 
N obtain estimates 

(5-2) \\U\\ H s Ut < C(\\P;U\\ HS; ^ m+1 + \\U\\ H - N ), 

with analogous results for s^. r . 

Now, as s^'*(a) = —s^(a) + (m — 1), if one chooses Sf tr as above, then 
one can take s* ast = — Sft r + (m — 1): with this choice, 

i.e. the space on the left hand side of (|5.1|) is dual to the (non-residual) 



space on the right hand side of (5.2), and (non-residual) the space on the 
right hand side of (5.1) is dual to the space on the left hand side of (5.2). 
Taking N sufficiently large such that the inclusions of the spaces on the left 
hand side of (5.1), resp. (5.2), into H~ N are compact, this implies Fredholm 
properties at once for P a and P*, with a slight change in the spaces as 
follows. Let 

ysftr-m+l _ j^Sft r -m+l ^ s ftr _ fJJ g ffS . p^jj g -ys ftr -m+lj 

(note that the last statement in the definition of X S{tI depends on the prin- 
cipal symbol of P a only, which is independent of a). 

Thus, we finally have the following, which follows from Propositions 2.2 
and 2.3 of [20J together with the propagation of singularities in variable order 
Sobolev spaces away from radial points (Proposition A.l in the appendix). 

Proposition 5.1. The family of maps P a enjoys the following properties: 
(1) 

P a : X SftT — > y Sftr_1 , P* : Af s p ast — > y s l^~ l 
are Fredholm. 

(2) If Pa depends holomorphically on a, P a is a holomorphic Fredholm 
family on these spaces with Sft r |A± = s ± ^ n 

(5.3) C S+)S _ = {a € C : s + < s + {a), s_ > s~(cr)}, 

while P* is antiholomorphic in the same region. 

(3) If P a is invertible (or if simply u G X Sftr , f G > ?Sftr_1 ; P a u = f), and 
WF(/) is disjoint from the radial set from which we are propagating 
the estimates, i.e., A~ , then WF(P~ 1 /) is also disjoint from this 
radial set. 

(4) If f is C°° , then WF(i- > ~ 1 /) is a subset of the radial set towards 
which we are propagating the estimates, i.e. A + . For the adjoint, 
corresponding to propagation in the opposite direction, we have 

WF((P CT *)- 1 /) C A- 
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when f is C°° . 

For the semiclassical problem, a natural assumption is non-trapping, i.e. 
all semiclassical bicharacteristics in £± apart from those in the radial sets, 
in Eft ± are required to tend to L + in the forward direction and L~ in the 
backward direction in £+, while the directions are reversed in £_. Here 
is the image of in S*X under the quotient map, and one considers 
S*X as the boundary of the radial compactification of the fibers of T*X. In 



particular, the non-trapping assumptions on M made in Section 3.2 imply 
that the operator P a is semiclassically non-trapping. 

Under this assumption, one has non-trapping semiclassical estimates (ana- 
logues of hyperbolic estimates, i.e. with a loss of h relative to elliptic es- 
timates), which, in the non-semi-classical language employed here, corre- 
sponds to an understanding of asymptotics as Recr -> oo. The following is 
proved in the same way as Theorem 2.15 of |20j. 

Proposition 5.2. If the non-trapping hypothesis holds, then: 

(1) P^ 1 has finitely many poles in each strip 

a < Imcr < b. 

(2) For all a, b, there exists C such that 



5 ^; 

ll^lyfe-i^ftr <C(Recr) _1 

on 

a < Imcr < b. 



Here the spaces with the subscript \o~\~ 1 denote semiclassical Sobolev 
spaces with the semiclassical parameter given by h = \a\~ ; we refer the 
reader to §2.8 of [20] for details. 

6. CONORMALITY OF COEFFICIENTS 

In this section we show that the coefficients in the asymptotic expansion 
which will appear in the sequel are in fact classical conormal distributions 
with a very explicit singular structure. 

For Y a connected component of {v = 0} in X, such as S+, let Aid 
denote the module of first order pseudodifferential operators on X = dM 
with principal symbol vanishing on N*Y. In particular, ^°(X) C Aid- Note 
that vector fields tangent to Y in fact lie in Mq, and indeed if A £ Aid, 
then because N*Y is locally defined by v = 0, r\ = 0, and 7 is elliptic on 
it, <ti(A) = aovy + Yl a jVji where dj G S°. In particular, A = Aq(vD v ) + 
Y^AjDy.+A', where Aj , A' G ^°(X), and so vD v ,D Vj , and Id generate Ma 
as a ^I /0 (X)-module. 

Belo w we work with the L 2 -based conormal spaces I^iX) defined m 



Section 2.3 above. Recall that u £ 7 (s) (X) means that u £ H S (X) and 
A\ . . . A k u G H S (X) for all k G N and Aj £ Mq. Thus is preserved by 
elements of Aig, while elements of ty k (X) map P s ' to l( s ~ k \ In particular, 
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when restricted to a product neighborhood of Y, elements of can be 
considered as C°° functions on Y with values in distributions on (—5, 5) 
which are conormal to {v = 0}, i.e., l( s \N*Y) = C°°{Y; j( s )(iV*{0})). 

We also recall the standard conormal spaces, defined using the L°°-based 
symbol spaces: a £ S k (Y x (—6, 6);M) if a is a compactly supported (in the 
(y, v) variables) and smooth (in all variables) and satisfies the estimates 



DyD e v D^ a 



k-N 



Elements of I r (N*Y) are then given as oscillatory integrals (which are es- 
sentially partial Fourier transforms), i.e., u £ I r (N*Y) if and only if 

u = J e w "<a(v,y,~f) d 7 with a £ S r+{n ~ 3)/4 , modulo C°°, 

i.e., a 6 S k corresponds to u 6 I k ~( n ~ 3 )/ 4: (N*Y). Since a £ S k corresponds 
to a lying in the weighted 1? space (^f) kJrl / 2+t L 2 for e > 0, 

jfe-(n-3)/4^y*y) c pj j(-k-l/2-e)^ N *Yy 

Note that N*Y corresponds to v = in this parameterization, and so the 
principal symbol is identified with an elliptic multiple of a| w= o- 

Now, if a is homogeneous outside a compact set in 7 (and a is independent 
of v near v = 0), one regards it for convenience as a homogeneous function 
on Y x (K\ {0}), and then a basis of such functions of degree k over C°°(Y) 
is given by 7 K times the characteristic function of (0,oo) 7 , resp. (— oo,0) 7 , 
which we denote 7^. If k is not a negative integer, one can go further, 
and consider the homogeneous distributions x±{l) 011 ^ ( or Y x K m our 
setting) defined by (the analytic continuation in n, from k > — 1, when they 
are locally L , of) 

= r(i + #c)' 

The inverse Fourier transform of these distributions are elliptic multiples of 

(see Section 7.1 of [9|); these are thus a basis for I^(N*Y)/I k f 1 (N*Y) for 
k = k - (n - 3)/4 over C°°(y). (The "cl" subscript stands for "classical" 
and refers to conormal distributions whose symbols have polyhomogeneous 
asymptotic expansions.) For negative integers k = —k, one must be more 
careful in describing a basis, as x± k is then supported at the origin. We 
instead simply consider directly the inverse Fourier transform of ^(±7) 7 _fc 
where ip is a smooth function equal to for 7 < 1 and 1 for 7 > 2. The result 
is a fc'th antiderivative of V , (=t 7 ), whose inverse Fourier transform differs by 
a smooth function from a multiple of (v ± «0) _1 ; hence differs by a smooth 
function from a multiple of 

(6.1) v ±lo K = (vizOr^log^iiO). 
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Note that these are no longer homogeneous distributions. (We also remark 
that when k is a negative integer of course we may also write more simply 

v ±lo K = v ~ l ~ K ( lo § M ± 
with H the Heaviside function; however it is more convenient to stick with 
the consistent notation offered by the expression as (6.1).) 
At first we consider a general operator 

(6.2) p = vD 2 v + aD v + Q, Q G Mg, a G C°°(Y); 

since D v is elliptic on N*Y, we in particular have P G A4g 1 ^ 1 (dM). Note 
in particular that the operator family P a defined by (3.6) has the form (6.2) 
by (3.9), hence the results here apply if P a u G C°°. 

Lemma 6.1. // Pu = f G I^(N*Y), u G l( s \N*Y), then 

(6.3) u = g + v 1 + i r+9-v 1 _^ a + u, 

with g± G C°°(Y) (pulled back via a local product decomposition) and u G 
l(*+i-e)(]V*F) for all e > 0. 

Remark 6.2. If t^-To" £ /( s+1 ~ e ) for all e > 0, then the conclusion is simply 
u G l( s+1 ~ e \N*Y). On the other hand, if u±7(T ^ t ' ien tne conclusion 
is g± = 0, and thus again u G l( s+1_e )(iV*y). 

If / G C°°(9M), iterative use of the lemma yields a full expansion of 
u, provided we replace g± by appropriate functions g± with P(g±v}^Q a ) G 
C°°(X) (see Lemma 6.4 below). 

Proof. We may assume that u is supported in a product neighborhood of Y, 
identified as (— 5, 5) v x Y, since if % £ C°°(<9M) is compactly supported in 
such a neighborhood and is identically 1 near Y, then WF' ([x, P])nN*Y = 
0, so [x, P]u G C°°(X) and thus P(xu) G C°°(X) as well. 

Note that vD 2 = D v vD v +iD v . Thus, if G G VP (Jf) is a parametrix for 
D v near iV*y (where D v is elliptic), applying G to yields 

+ {% + a) + GQ)u G I (s+1) (iV*y). 

Since u G /W, Q-u G I^- 1 and thus GQu G j( s+1 ), so we have 

(vD v + i + a)u£ = C°°(Y; I^ S+1 \N* {0})) . 

With / a compact interval, let I^\n*{0}) denote the sum of elements of 
/W(_/V*{0}) supported in / and Schwartz functions on R. Now, note that the 
Fourier transform on R maps elements of I^\n*{0}) to L 2 -based symbols. 
More precisely, if S^' is the set of smooth functions (j) on R 7 such that 

( 7 D y f<p G L 2,e = < 7 >^L 2 
for all N G N, then the Fourier transform is an isomorphism Ig\N*{0}) — > 
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Taking the partial Fourier transform, J 7 , in the interval variable, v, yields 

(-7-D 7 + 1% + a)Fu = (-Ztyy + i + a)Fu G C°°(Y; S {s+1) ). 

Now, to analyze the behavior of Tu at infinity, we conjugate the differential 
operator by 7 _2+iQ on K \ {0}, where 

7 2— (_ 7 £, 7 + 2i + a ) 7 - 2+4Q = - 7 D 7J 

so one has 

-jD^-^Fu) = 7 2 - iQ (- 7 D 7 + 2i + a)Ju G C°°(Y; ^- 1 " Im °)) ) 
and thus 

(6.4) J D 7 ( 7 2 -^n) GC°°(y;5 (s - Ima) [l,oo)). 

Note that due to the presence of e > in the statement of the lemma, 
we may assume that s — Ima 7^ 1/2; this simplifies some formulae below 
(otherwise one would have logarithmic terms). 

Now, if b G 5^([1, 00)), £ < 1/2, then the indefinite integral of b given by 

/7 
6(77) tfy, 

satisfies (by Cauchy-Schwarz) 

k 7 )i< (j"\rj\ 2e \b( V )\ 2 dS /2 ( r\v\- 2i dj 1/2 

<c||% 2 ,, (i + | 7 |^ 

Thus c G L 2,f_1 ~ e for all e > 0, and as L> 7 c = 6, c G S^ 1-6 ). (Note that 
constants are in S^ -1- ^ since £ < 1/2.) 

Returning now to u described by (6.4) above and setting £ = s — Ima, 
we see that 

1 2 ~^pu = w G S( s - Im «-!^). 

On the other hand, if i = s — Ima > 1/2, then S^' C L 1 , and if we define 
the indefinite integral as 

poo 

c(7) = - / Kv) drj, 

then, by Cauchy-Schwarz, 

/ roo \ 1/2 / roo \ 1/2 

|c( 7 )|<(J |r?| 2 ^W| 2 ^J (jf M- M A/J <C||6|| L2 ,| 7 |"^, 
so c G S 1 ^ -1-6 '. Then, writing 



00 



we deduce that 



<f-**Fu = 5+ + 5+ G C°°(Y), w G S'(«-Ima-l- £ ) ! 
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and thus 

Fu\ 1>x = g+r 2+ia + w+, w + G 
A similar calculation applies to J-"n| 7< _i, yielding 

JT U | 7< _ 1 = ff _(_ 7 )-2-Ha + tw _ ) ^_ e 5-(«+l-e). 

In summary, if ip + is supported in (l,oo), identically 1 on [2, oo), and 
V'-(t) = then 

fu = g + i; +1 - 2+ia + 5 _V^(-7)" 2+Ja + w, w G S (s+1 " e) . 

Now, the inverse partial Fourier transform of to is in /( s + 1_e ), so it remains 
to deal with the other terms. Changing these by a compactly supported 
distribution does not affect their singularities, so we can replace these by 
the homogeneous distributions 7^ 2+lct for a local description of the inverse 
partial Fourier transform. But the inverse Fourier transforms of the latter 
are v^ a , so we conclude that 

as claimed. □ 



Although the following corollary follows directly from [20j , we give a proof 
using Lemma |6.1[ 



Corollary 6.3. If Pu = f G C°°{X), u G l( s °\N*Y), s > 3/2+Ima, 
then u G C°°(X). 

Proof. Let s = sup{s : u G I^(N*Y)}, so so > 3/2 + Ima (possibly 
so = +oo); if so = +oo, then we are done as PlseR-^ = C°°(X). Thus, 
u G j(So-e) (AT*y) f or a n e > 0. By Lemma 



6.1 



u = g + v 1 + ™ + g„v 1 _™ + u, 



with g± G C°°(y) (pulled back via a local product decomposition) and u G 
j(s +i-e)( N * Y } f or a u e > o. For all e > 0, u G /( So+1 - e )(A^*y), which is a 
subset of /( 3 / 2 + ImQ )(AT*y) for sufficiently small e > 0. On the other hand the 
sum of the first two terms is not in /( 3 / 2 + Ima )(Af*Y) unless g± vanish. Since 
u G j(3/2+im a )^ N *Y^ g± must vanish, and thus u = u G j( if o+i-e)(JV*Y') 
for all e > 0, contradicting the definition of so- Thus, s~o = +oo, completing 
the proof. □ 

Next, under the assumption that a is constant, we show that distribu- 
tions such as those in the first two terms on the right hand side of the 



equation (6.3) can be modified to elements of the nullspace of P modulo 
C°°(X). 

Lemma 6.4. Suppose a G C is a constant, 1 — la is not an integer, and 
g G C°°(Y). Then there exist u± = g±v^ a G n e>0 J (3/2+ImQ! " e) (-^*> r ) ) with 
g± G C°°(X) such that g±\y = g and Pu± G C°°(X). 
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Remark 6.5. If 1 — tot. is an integer, the proof below still proves a slightly 
different result: logarithmic terms appear. Indeed, if 1 — ta is a nonnegative 
integer, then logarithmic terms appear from the definition of v^q*. If it is 
a negative integer, say, 1 — la = —r < —1, then an additional logarithmic 
term is incurred at the r-th step in the expansion. Observe that in this case, 
the additional log term appears at the step in the expansion corresponding 
to v± l0 for nonnegative k. 

It is more straightforward to state it as follows: u± is a classical conormal 
distribution of the appropriate order, with principal symbol the same as that 

Remark 6.6. A similar expansion can be obtained in general, without as- 
suming that a is a constant. This is similar to the treatment of generalized 
Coulomb type spherical waves in |18j . 

Proof. We suppose first that 1 — la is not an integer. As the indicial roots 
associated to the ordinary differential operator vD% + aD v are and 1 — la, 
for h G C°°(X), 

Pv k v l ±™h = v k vg%w, w G C°°(X),u;|y = c(k)h\ Y , 

with c(0) = 0, c{k) ^ for k ^ 0. (We suppress the dependence of c(k) on 
a.) Correspondingly, given g, consider first h± t o G C°°(X) with ho\y = 9- 
Then 

with w\y = 0, so in fact 

Now, in general, for k ^ 0, given w±^ G C°°(X), one can let h±^ = 
—c(k) w± i k, and then 

Pv k v^ a h± jk + v k v±™w± t k = v k v^w± tk 

with w± fe|y = 0, thus the right hand side is of the form v k+l (v±i0)~ ia w±^+i- 
Correspondingly, one can proceed inductively and construct h± & with 

Pv^ a h± lk = v k+1 vg^w ±>k+1 , 

with w±„k+i G C°°(X), namely h±^ = Ylij=o v ^±j works. More generally, 
one can asymptotically sum the series X]j=o v3 h±,ji i- e -> construct a function 
h± which differs from Ylj=o v ^±,j a C°° function vanishing to order 
then 

Pv l ^ a h ± = v k+1 v^W ± , k+1 

for every k for some W± k+i £ C°°(X), thus the right hand side is C°°, 
completing the proof. 

If 1 — la is a non-negative integer, then the iterative construction requires 
including another logarithmic term owing to the logarithmic term in v± ia . 
If 1— ia = r < —1 is a negative integer, then the iterative construction breaks 
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down when finding the coefficient of v r v±Jf*, as in this setting c(r) = 0. The 

proof goes through nearly as stated once we also include terms of the form 
v r+k v v-w log ( v ± l0 ) for k > o □ 



In addition to knowing that we may formally parametrize elements in 
the approximate nullspace by functions on Y, we will need to know how 
to formally solve certain inhomogeneous equations with specified conormal 
right-hand sides. For the following lemma, we assume that Q is a differ- 
ential operator in the module .M|, although it holds (with a slightly more 
complicated proof) if Q is pseudodifferential. Note that for our operator P a , 
Q is in fact differential. 

Lemma 6.7. Suppose a € C is constant and Q is a differential operator 
in Mq. Let h £ C°°{X) and let m be a nonnegative integer. If ia is not a 
strictly positive integer, then there exist g±, . . . g± +1 6 C°°(X) such that the 
functions 

m+1 
m'=0 

solve 

p u± = hv£g Iog(« ± ^0) m + u 

with u G C°°. If 1 — ia = —k < is a non-positive integer, then the same 
statement is true with u± replaced by u± + w±, where w± has the following 
form with g± +2 a smooth function 

w± = fff 2 /- 1 ^" log(u ± i0) m+2 . 



If 1 — ia is a non-negative integer, then there is an additional log term 
implicit in the formula owing to the definition of v^q° . 

Proof. The proof is similar to that of Lemma |6.4| Indeed, as the indicial 
roots of vD"l + aD v are and 1 — ia, for g £ C°°(X), 

max(m',2) 

(6.5) P (v^^ilogiv ± iV)T'g) = vkv Zlo(}og(v ± *0)) ra '-^W, 

w (t) eC °°(X), w^\Y = c{k,m',£)g\ Y . 

Note that 

Here c(0, m',0) = 0. If ia is not a positive integer, then all other 
c(k,m',£) 7^ 0. If ia = 1, then c(0,m',l) = as well and the other 
c(k,m',£) 7^ 0. If ia = r is a positive integer, then c(r — l,m',0) = 
with all other c(k,m',£) ^ 0. 
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We start by assuming that ta is not a positive integer. Consider first 
g m+l g C oo( X y We have 

max(m+l,2) 

= E «°*£o (M« ± *o)r +1 -^V+i< 

£=0 

with u4\ m+1 = c^m+M)^ 1 . Asc(0,m+1,0) = 0, and e(0,m+l, 1) / 
0, we set c^ 1 = c(0, m + 1, 
Given g±' Q G C°°, we have 

Pt;°4X(log(«±,0)) m '^; 

max(m',2) 

e=o 

For < m' < m, we then set 

, __c(0X + 2 1 2) (2) 
5± '°" c(0,m+l,l) ± -»' m '+ 2 ' 

Applying P, all terms other than those involving tOj? m' cance l a ^ All 
^+ o m' van i sn a * ^ an d so in fact 
m+l 



P E ^ 4;o°(log(^ ± ^)) m '< - v££ Qog(v ± tO))™* 
m'=0 
m+l 



m'=0 

Now, in general, for / 0, given w±\ m i S C°°(A) for < m' < m + 1, 
we set 

0±f* = ~ c (jfc m / )0 ) (™±.M* # + C ( fc ' m ' + ^iim'+l + C ( fc > m ' + 2 ' 2 ) W ±i,m'+ 2 ) » 



where uA . , are the coefficients in equation (6.5) with applied to g = g™ . 
Applying P, we see 

m+l m+l 

P E v k v^(log{v±zO)r'g^ k + E «^£o 

m'=0 m'=0 

-^(logCuitO))^ 

m+l 



•"±10 w ±,k,m'> 
m'=0 
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where w±^ k ,\y = 0, and so the right hand side is of the form 

m+l 



53 V k+1 V±lo™±,k+l,mi- 
m'=0 

We can thus proceed inductively and construct g±' k with 



m+l 



p 53 4X(iog(« ± M)r'~g±' k - ^(Mv ± ri))^ 

m'=0 

m+l 

53 u fc+1 u3S'is± I fc + i ir . 



m+l 

"/ +1 ^0*±,fc+l,m', 
m'=0 

with ^ ±ife+ljm . G C°°(X). (Namely, ^ = £* =0 ^S^ works.) 

We now asymptotically sum the series Yl'jLo v " ' 9±,j *° construct a function 
g± differing from each ^j=o ^ 9±j by a smooth function vanishing to order 
k + 1, and then 
m+l 

P 53 t43°(log(« ± *0)r'<#' - «£?(log(i> ± iO))™* g C°°(X), 

m'=0 

completing the proof. 

If 1 — la = — k is a non-positive integer, the iteration proceeds nearly as 
before, but at the expense of an additional log term at the k-th. coefficient. 
(For example, if k = 0, then c(0,m', 1) also vanishes and so an additional 
log term is needed to find the first coefficient.) □ 



We now combine Lemma 6.1 and Lemma 6.4 to obtain a complete as 



ymptotic expansion of elements of the nullspace of P modulo C°°(X). 

Proposition 6.8. Suppose a G C is a constant. IfuG I^ S \N*Y) for some 
s and Pu G C°°(X), then there exist g± G C°°(X) and u G C°°(X) such that 

u = 9+vlT + + u. 

See Remark \6.h^ if 1 — ta is an integer. 

Proof. Let so = sup{s : u G 1^ (N*Y)} (the set on the right is non-empty 
by hypothesis), so either so = +oo, and then HseM ^ = C°°(X) shows that 
the conclusion holds with g± = 0, or sq G R is finite, and then u G j( So— e ) 



for all e > 0. By Lemma 6.1, there exist g± G C°°(Y) so that 

u=~9+v 1 ^ a +~9-V 1 ~v a +u, 

with v! G l( so + 1-<5 ) (]V* Y) for all 5 > 0. Here the first two terms are in 

r\6>o l{3/2+lma ~ 5) ( N * Y ) but not in I^I 2+lraa \N*Y) unless g ± vanish; by 
the assumption on sq, 3/2 + Ima = sq and g± cannot both vani sh. Let 



g± G C°°{X), u ±ef]s>o I^/ 2+lma ~ s \N*Y) be given by Lemma 6.4 with g± 
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in place of g. Thus, Pu± G C°°(X), hence P(u — u + 
Further, 



u-) G C°°(X). 



u — u + 



U- 



and 



pi j(5/2+Ima-5)^jy*y^ 

Thus, u - u+ - u_ G f\>o j( 5 / 2 + In » 0[ -*)(JV*l r ). By Corollary 
G C ,00 (X), completing the proof of the proposition. 



6.3 



u — u + 



□ 



In our setting, where by equation (3.8) 



a = a 



this gives: 

Corollary 6.9. If u € I^°°\N*Y) and P a u G C°°{X), then there exist 
g±€C°°(X) and u£C°°(X) such that 

(6.6) U = g + vZg + g-Vl% + u. 

Again, see Remark\6.5\if —la is an integer. 



Note that u as in the corollary lies in £f 1 / 2 + Imo "- e f or a \\ e > 0, but not 
in ffi/s+imo- un j ess g±\y vanish. Thus, for s and a corresponding to the 
region (5.3), this lies in H s , the target space of {Pa)Ytv > as expected - and 
this containment is sharp in so far as it would fail whenever g±\y do not 
vanish if the inequality in (5.3) is replaced by equality. 

Finally, we now use Lemma 6.7 to deduce the structure of solutions to 
certain inhomogeneous equations with conormal right hand side: 

Proposition 6.10. If u G I^°°\N*Y) and 

P a u£vZ%- 1 log(v±i0) m C oo (X) J 

then there exist g±' G C°°(X) (for ml = 0, . . . m + 1) and u G C°°(X) such 
that 



(6.7) 



u 



m+l 
m'=0 



See Remark\6.5\ and Proposition 6.1 if —io~ is an integer. 



Proof. By Lemma 6.7, we may find a function w of the form (6.7) so that 

P a w- P a ueC°°(X), 

with the leading term having the claimed form. As the function w is also 
conormal, w — u is conormal, and so we may apply Corollary 6.9 to finish 
the proof. □ 
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7. The connection between P a and asymptotically hyperbolic 

and de Sitter spaces 

In this section we identify the action of P" 1 on the caps C± and in 
the equatorial region Co as in [201 Sections 3.3 and 4]. Recall that P~ l 
propagates regularity from S- to 5+; this means that the behavior at C_ 
and Co is what is studied in detail in the above references, with the behavior 
of P~ l at Co and C + corresponding to the adjoint operator in that paper. 

So we consider the operators 

L ff ,± = -A fe± + (a 2 + (n ~ 2 H + vX{±a) + vV(±a), 



from (3.10) (note the sign switch in a relative to (3.10) to keep the behavior 
for -Lo-,+ and L a - similar in terms of Im a > being the physical half-plane), 
with k± asymptotically hyperbolic metrics, V a smooth potential and X a 
vector field tangent to v = 0, as well as 

L afi = n feo + (a 2 + (n ~ 2 H + vX[a) + vV(a) 



from (3.11 ), with V, X as above (\v\ = —v being a defining function for Co). 

Since L CTj o is an asymptotically de Sitter operator as in |22| . it has a for- 
ward solution operator 1Zc (cr) propagating towards S+, i.e. if / 6 C£°(Co), 
u = TZc {cr)f is the unique solution of L a pu = f with Lq vanishing near S-. 
On the other hand, L a ± are non-self- adjoint perturbations of the asymp- 
totically hyperbolic operator — Afc ± + (a 2 + ^"~ 2 - > ), as in [13], with the 
perturbation being non-trapping in the high energy sense. In particular, 
L a ± : Hq[C±) — > L 2 (C±), Imo" > 0, is an analytic Fredholm family. Since 
in general we do not have automatic invertibility for such perturbations 
without appeal to the large parameter behavior, which is only understood 
from the perspective of the extended operator, we need a proposition. Note 
that if vX(a) + vV(a) vanishes then the invertibility of L a is automatic 
when Imo" > 0, Imu 2 7^ as A^ ± is self-adjoint. 

We then have the following lemma: 

Lemma 7.1. Suppose that P a : X SftI ->■ y s ^- 1 and P a : X s ^ st ->■ 37 s p^t-i 
are invertible for a certain a G C with Im a > 0. Then L a> ± : Hq(C±) —> 
L 2 {C±) is invertible. 

Remark 7.2. While we handle the invertibility within our framework, an 
alternative would be the complex absorption framework used in [p5]; the 
absorption would be placed in v < — e for some e > 0. 

Proof. As already remarked, L CT) ± : Hq(C±) — > L 2 (C±) is Fredholm, so we 
only need to show that Ker L a< ± and Ker L* _j_ are trivial. Using [13] , first any 
element of Ker L CTi ± is in Hq°(C±) by elliptic regularity in the 0-calculus, and 
indeed using the parametrix construction, they are in v ™ 2 + 4 ~5~C°°(C±), 



ASYMPTOTICS OF RADIATION FIELDS 



11 



1 n 1(7 Q 

while any element of KerL* ± is in v~? + ± ~~2~C°°(C±). In particular, for 







L a — , for any element U- of the kernel, we can extend v% 4 2 «_ to an 



element u otC°°{X). Then / = P a u is supported in C + U Co by (3.10), hence 
P" 1 / is also supported in C+ U Co, so u = u — P" 1 / solves P a u = and 
u\c~ = u -- Since Ker P a is trivial by assumption, u, and thus U-, vanish. A 
similar argument applies to elements of Ker L* , as P* : Af s p ast — )• 3^ s p ast_1 is 

also invertible; in that case for an element ii+ of the kernel vz~i + ~u+ to 
an element u of C°°(A A ) and apply (P*) _1 to the result. Finally, for Ker L a + 
and Ker L* _ we switch the direction of propagation for the inverse P^ 1 , i.e. 
we consider 

and then completely analogous arguments apply as the roles of C+ and C_ 
are simply reversed. □ 

We now make the connection between P~ l and the operators on the C± 
and C . K c _(.) is regular at -a , lma < 0. If / G C C °°(C_) C C°°(X), 



then P 1 / is smooth on C_ by Proposition 5.1 By (13.10), on C 



U = U 2 + 4+ 2 (P ffo i (t;4+ 2 + 2/))| c _ 
1 1 n 1 'Q-Q 



solves = /, and u G v~^ + ^ + ^rC°°(C^) C L§(C_) (with L§(C_) 

being the asymptotically hyperbolic L 2 space), with the inclusion holding 
as Re(icio) = — Imcro > 0. Thus, 



P£( V $+*?+$f))\ c _=Kc-{-<T )f, 



1 1 « 1 JCT / „_1 , n 1 
2" r 4" r 2 



since 7£c_ ( — Co)/ is the unique L 2 (relative to the asymptotically hyperbolic 
metric) solution of L_ CTOi _u = f. By the meromorphy of both sides, the 
formula is then valid at all <jq (regardless of the sign of Imo"o) at which 
P~ v is regular. Indeed, by the same argument, for / G C°°(C_) such that 
t;4+^+2 f £ C°°(C_), and thus has an extension / = E(vi ^ zf) to an 
element of C°°(X), one still has 

U -2+4+-(P-^( U 4+-+2/))|c_ =n C .{-(T )f. 

Further, if / G C£°(CoUC + ) and P" 1 is regular at 00, then P~^ f vanishes 
in C_. Indeed, first if Imcro < 0, and if it did not vanish, then PjJ, 1 / would 
yield an I? solution of L-^—u = in C_ and correspond to a pole of 
7^C-(— co) and hence to a pole of P" 1 by the above discussion. As P&} f is 



smooth in a neighborhood of S- by Proposition 5.1 P^, 1 / must vanish to 
infinite order at S-. Turning to the region Co, the Carleman-type estimates 
in [22l Proposition 5.3] (see also [201 Section 4]) imply that P^ 1 / must 
vanish in a neighborhood of S-. In particular, P~^- (applied to such /, with 
the result restricted to Cq) must be a conjugate of the forward fundamental 



solution of the operator in equation (3.11 ) (applied to f\c ), which we denote 
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lZc (cr). Indeed, again a simple generalization shows the same conclusion 
when one merely has / 6 C°°(Cq U C+), with the dot denoting infinite order 
vanishing at the boundary of this set, namely S-. 

Finally, if / € C£°(C + ), or indeed / G C°°{C + ) then the above discussion 
implies that P^f vanishes in C_ and Cq. Moreover, the expansion of 



i-— — 

\v 2 4 " 



Corollary 6.9 implies that in fact P ao is a conjugate of TZc+(o~)- 

In particular, the above discussion proves the following proposition: 

Proposition 7.3. If P^ 1 is regular at a, then it has the following "block 
structure" (here the rows and columns correspond to support in C+, Cq, and 
C-): 

^■c+(o")hl 242 * * 

. 1 to %a_ 1 i to I 2^ 

|l>| 2 4 2 T^-C^o V ") I V| 2 4 2 * 

. . 1 n icr 1 I TO I icr 

o o |u| 5- 4-^7e c _ (— cr)|i;| 2 +4 

in the strong sense that if P~ l is applied to a C°° function on X , the re- 
striction of the result to C_ is given by the lower right block, 

. 1 n %a_ 1 ■ rt | i£ 

1 17 1 2 4 2 7£f7_(— C)M 2 4 2 ) 

if P^ 1 is applied to a C°° function supported in Cq U C+, t/ie restriction of 

I - 1 ti io_ . ,1 I n l Iff / J 7 

me result to Cq is given by \v \ 2 4 2 7cc ((t)|u| 2^4^ 2 (ana ims result van- 
ishes in C_, w/iiZe finally if P^ 1 is applied to a C°° function supported in C+, 
me restriction of the result to C+ is given by \v\ 2 * 2 TZc + {&)\v\ 2 * 2 
(and the result vanishes in C_ U Cq). 

By our non-trapping assumption on the null- geodesies of g, — A& ± + 
vX(a) + uV(<r) + cr 2 + (n — 2) 2 /4 is a semi-classically non-trapping operator 
and thus the following proposition (which follows from, e.g., the work of 
Vasy [20J) holds. 

Proposition 7.4 (cf. |1Q|, Theorem 4.7]). Consider the operators 
L a>± = -A k± + (a 2 + (n ~ 2 H + vX(±a) + vV(±a), 



with k± asymptotically hyperbolic metrics, V a smooth potential and X a 
vector field tangent to v = 0. If L a ± is semiclassically non-trapping, then it 
has a meromorphic inverse TZc ± (cr) with finite rank poles, is holomorphic for 
Imo" S> 0, and has only finitely many poles in each strip C\ < Im a < Ci. 
Moreover, non-trapping estimates hold in each strip Im a > —C for large 
Re a: 

\\TZ c± (a)f\\ HS <C \\f\\ H s-i 

Moreover, if L a ± has no I? "eigenvalues" (with respect to the metric 
k±), i.e. no element of L 2 lies in Ker L a ,±, then all poles uq ofIZc±(o~) have 
Imffo < 0. 
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Proof. The bounded a properties were already explained above. The high 
energy estimates then follow from those for P„ l . Since P a has index zero 
as a Fredholm problem, its invertibility amounts to having a trivial kernel. 
Since an element of KerP^ restricts to a C°° function on C_ U Co, thus 
«i _ 2+ ! 2 times the restriction to C_ is an element of ui - 2 + *2C°°(C_), in 
view of the asymptotically hyperbolic metric on C_ this gives an element 
of L 2 if Im a < 0, so under the assumption of no I? "eigenvalues" all poles 
o"o of TZq_ (c) indeed have Im o~o < 0. On the other hand, by Sectionj6j an 
element of KerP CT whose support is disjoint from C_ is supported in C+, and 
restricted to C+ it has the asymptotic form v~ ta C°°(C+), and thus i;4~2+ l 2 
times the restriction to C+ is an element of v *~ 2^*2 C°°(C_), and thus is in 
L 2 if Im a > 0, under the assumption of no L 2 "eigenvalues" all poles o~o of 
T^c+ ip) indeed have Im ctq > 0. □ 



Remark 7.5. Proposition 7.3 implies that the poles of TZc + (c) yield poles 
of P^ 1 - A partial converse is true as well. If do is a pole of P~ l so that 
the corresponding resonant dual state has support intersecting X \ C_ , then 
either oq is a pole of 7lc+ (c) or the corresponding resonant state is supported 
at S + (see [201 Remark 4.6] for more details). Such poles may occur only 
for o"o a pure imaginary negative integer. In other words, the relevant poles 
of P~ l are either poles of !Zc+ (o*) or have state supported at S+ (and hence 
are differentiated delta functions in v). We remark that such states occur 
in even-dimensional Minkowski space, where — % is a pole of P~ l in 2- and 
4-dimensions. 



8. Structure of the poles of P a 1 

While the results in the previous section fully address the structure of 
nullspace of P a , knowledge of nullspace alone is clearly not sufficient to deal 
with the structure of the poles of P^ 1 - Even for a spectral family of the form 

(Po-^Id)- 1 , 



with Pq as in (6.2 ), the poles may of course be multiple owing to generalized 
eigenspaces; thus knowing that the nullspace of Po has a particular conormal 
form would in general permit the range of the polar part of the resolvent 
to have log terms. Here the situation is further complicated by the fact 
that our family P a is not of the form Pq — aid but rather has nontrivial 
dependence on a, so that we cannot even employ the usual machinery of 
Jordan decomposition. A careful analysis of the log terms will, however, be 
essential in order to see that excess log terms in our asymptotic expansion 



(1.1) do not spoil the restriction of the rescaled solution to the front face 



of the radiation field blowup, which we know a priori must be smooth (cf. 



^3.7). In this section we demonstrate (among other things) that the top- 



order terms with log p are balanced by terms containing log v in such a way 
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as to permit the solution to be smooth across the front facep*| We should 
emphasize that these log terms are typically not vanishing, and are still 
a relevant part of the expansion away from the interior of the front face. 
In particular, we prove the following proposition, which is an extension of 



Corollary 6.9 



Proposition 8.1. Let o~o be a pole of order k of the operator family 

p-i . ysttr-i x s a ^ 

and let 

(a - a )' k A k + {a- <7 Q + . . . + (a - o- Q )~ l A l + A 

denote the Laurent expansion near <7o, with Aq (locally) holomorphic. If 
f £ ^ftr- 1 vanishes in a neighborhood of C_, there are smooth functions 
(j)± t i, . . . <ft±,k so that for < I < k — 1, Aef has an asymptotic expansion^\ 



Ak-ef 



llO 



+ v~_ 



^i-^(log(« + i 0))^ +ifc _ ( ,_ i) 



3=0 



i0 



^L^.0 O g(«-,o)^_^_ i) 



3=0 



J- 



+ v 



-UTQ + l 



(log v) 



If—icro is a non-negative integer, then there are smooth functions (f>i, . . . , <j>k 
so that Ak-ef has a similar expansion in terms of the distributions v^_ tcr ° = 

A k -ef = < CT0 Yl ^f( lQ g M)Vk-(*-;) + 0(v-^ +1 (log \v\) e ). 

3=0 h 

Remark 8.2. This proposition serves two purposes. The first is to show that 
Laurent coefficients have asymptotic expansions at v = 0, while the second 
is to show that the leading terms in this expansion have a specific form. 
This form is later used to show that the terms of the form log p cancel at 
the radiation field face so that the radiation field may be defined. 

The additional logarithmic terms occurring at imaginary integers in Propo- 



sition 6.7 would in general disrupt the form of this expansion, but we use 
the support of the states to conclude that in fact it has the desired form. 

One could also write the entire expansion in terms of H(v)v~ ta ° even 
if — igq is not a positive integer. To do this, we would have to include 
derivatives of delta functions if —ictq is a negative integer. (These occur 
even in the case of even dimensional Minkowski space.) 



1 *-$ 

In light of the smoothness of the solution across the front face, we expect all such 
log terms to be balanced in this manner, but we consider only the top-order terms (the 
terms affecting the radiation field) here. We have not undertaken this here, however. 

14 Although we use the notation 0(w 7 (logf) K ) here, the term in fact has a polyhomo- 
geneous expansion with index sets shifted from the "base" ones. 
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The proof requires the following lemma: 



Lemma 8.3. If f vanishes in a neighborhood of C_ then Ak-ef is supported 

inC + /or £ = 0,1,..., A; -10 

Proof. Near a pole cr$ of P^ 1 , we may write 

(8.1) P a = P + (a-a )P 1 + (a-a ) 2 P 2 , 

where Pq = Pr > Pi = D v + E, and P2 is a smooth function. Here E G 
A4q is a first order differential operator characteristic on N*S+. The proof 
relies on the following relationships between Pi and Aj , which holds because 

P.Pa 1 = I: 

(8.2) Po^fc = 

PiAfc + P 4fe-i = 

P 2 ^ fc _i + PxA k _ {i _ x) + P A Mi _ 2) = 0, i = 0, . . . , k - 3 

We first observe that A k _if vanishes near C_. Indeed, for £ = this 
follows from the Cauchy integral formula applied to (a — ao) k P~ 1 f and 



Proposition 7.3, while for £ > 0, it follows inductively from Proposition 7.3 



the Cauchy integral formula applied to 

£-1 

(a - a Q ) k+t p- l f - J> - a YA k _ 3 f. 
3=0 

To observe that A k _gf vanishes in Co, we again proceed inductively. For 
i = 0, as PoAkf = 0, Proposition 7.3 implies that it vanishes in a neighbor- 
hood of S- and hence in all of Co- If I > 0, the relationship (8.2) implies 
that P^Ak—if vanishes in CqUC_ and so A k -ef also vanishes in CqUC_. □ 



Indeed, 



Proof of Proposition 8. 1 . We rely on the structure of P a near S. 
recall from above that P a = D v (vD v + a) + Q, where Q G 
notation of Section [6]) is a differential operator. We rely on the form 
P(j near a pole o"o of P" 1 as well as the relationships (8.2) between Pj and 



Mj (in the 



1) of 



We start by assuming that 
on £. As / vanishes near C 
in C+, while Proposition 



5.1 



^o"o is not an i nteger and proceed by induction 
implies that A k f is supported 



8.3 



Lemma 

or, indeed, elliptic regularity) implies it is 



smooth away from the radial set A . We may thus apply a theorem of 
Haber-Vasy P Theorem 6.3] to conclude that in fact A k f £ I^-°°\A + ) = 



l(-°°)(N*S + ). In particular, Corollary 
functions <j>± k and ip so that 



6.9 



implies that there are smooth 



A k f 



J +i0 



15n 



This lemma implies that there are two types of "resonant" states. If the state is given 
by •}, then either <j) is supported in C+ or ip is supported in C_. See [201 Section 4.9], 
especially Remark 4.6, for more details. 
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By Lemma 8.3, A k f is supported in C+, so tj; vanishes to infinite order at 



S+ and may be absorbed into the other terms, i.e., 

Now suppose that the statement is true for < £' < I — 1. As PoA k _ef = 
-P 1 A k ^ e+1 f - P 2 A k _ e+2 f, we have 

PoA k „ e f = e E u ±r _1 ( io g( u ± w^r- [(7o^±,fc-«-i-i) + ^±,k-(/-2-i)] 

+ O(«- l,70 (logt;) £ ), 

where the O(t> _iCT0 (logu)^) in fact has an asymptotic expansion of a similar 
form. Observe that the right hand side is an element of l(~°°)(N*S+ ), so 
again Haber-Vasy implies that A k _if G l(~°°\N*S+). Proposition 
then implies that has a similar expansion, say 

l 

Ak-tf = E E v ±*o° ( l °s(v ± tO^Oij 
± j=0 

To determine the leading coefficients in the expansion, we calculate 
e-l 

P A k . e f = - EE ^°~ 1 ( 1 °g( u ± »0)) i C7" + 1) [(-^ )a ±J+1 + (j + 2)a ±)i+2 ] 
+ O(*r lCT0 (log^), 

where again the last term has an expansion. We now simply equate coeffi- 
cients, starting with the largest one. If j = £ — 1, we must have 

uro£a±j = — ^y Cro(j)±, k , 

(—i) e 

i.e., a± t i = ^jf-4>±, k - Now for j < £ — 1, we have 



ta-oO" + l)o±j+l - (i + 2)(i + l) 

(fO0±,jfc-(f-l-j) + 0±,fc— 2— j)) > 



i.e., a±j+i = /(? + 1)!^± fc-(£-i-j)- This determines a±, 2 , • • • ,a±,£, 



while a± i are given by Corollary 6.9 and are denoted <f>± k _£. 



We now consider when — ioq is an integer, in which case additional log- 
arithmic terms appear in Proposition |6. 10 If — ictq < 0, these additional 



logarithms are not in the leading order terms and so the results above still 
hold. For — i<jq > an integer, however, we must be a bit more careful and 



rely on Lemma 8.3 
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Let us assume for now that —ictq ^ 0. Indeed, we again proceed induc- 
tively. Consider first Akf- The same arguments as above imply that A k f 
has an expansion of the form 

A k f = v- + \l^ + + vZ7^- + ^ 

As Akf is supported in C+ and v±?o° = v~ lao \og(v ± tO) in this case, given 
the behavior of 0_ and 4> at 5+ is determined by the support condition. 
Indeed, we must have that 0_ = — <p + and <j> = —2iriv~ la °(f)+. In other 
words, there is a smooth function so that 

A k f = v-^H(v)^ k . 

Now suppose that the statement holds for Ak-pf for <£'<£ — 1. We 
then have that P$Ak— if must satisfy 

Po^W = - ^T^'-^Xlog M)'^- (a o 0*-(^i-i) + fa-V-*-])) 

3=0 J ' 

+ O(v- W0 (]og\v\) e ), 
where again the 0(w~ tCT °(log |v|) ) term has an expansion of a similar form. 



The Theorem of Haber-Vasy and Lemma 8.3 then imply that Ak-ef has an 
expansion 

i 

A k _J = Y,v~ la °(log\v\ya r 

3=0 

Applying Pq and equating coefficients finishes the proof in this case. 

Finally, if — ioq = 0, the same argument as in the case of — ioq > still 
works, but differentiating the j = term yields a 5{y) term. This term is 
no problem, as we still simply solve for its coefficient. This process yields 
an identical result. □ 



9. An asymptotic expansion 

In this section we detail the iteration scheme required to obtain a pre- 
liminary asymptotic expansion for (smooth) solutions w of O g w G C£°(M°) 
that vanish in a neighborhood of C_. Recall that in the notation above, 

□ 6 = p- 2 Ug. 

We start with a tempered solution w of Ow = f G C^°(M°) vanishing in a 
neighborhood of C_. We may replace w by x w i where x ' ls a cut-off function 
supported near the boundary. The right hand side is then still smooth and 
compactly supported, w still vanishes in a neighborhood of C_, and w is 
supported in a collar neighborhood of the boundary. As before, write 

L = p -(n-2)/2 p -2 np (n-2)/2 £ Diff 2 (M)> 

so that setting 
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we have 

Lu = g£ C C °°(M°). 

Now let N(L) denote the normal operator of L and set E = L — N(L); E 
thus measures the failure of L to be dilation-invariant in p. Thus, 

E 6 pDiffg(M). 

Because w is tempered, we know that w 6 p 7 -f/^°(M) for some sq and 



7. After possibly decreasing so so that 7 + so < 1/2, Corollary 4.6 implies 
that u> in fact has conormal regularity with respect to N*S+ relative to 
p' y H^ (M). By the form of G~ l given in equation (3.3), we note that the 



coefficient of D 2 in E is of the form 0{p 2 ) + O(pv), hence Ew G p^^Ht ^ . 

At this juncture, we remark on the mapping properties of the Mellin 
conjugate of E. To begin, we let R a be the family of operators satisfying 

M o E = R a o M; 

thus R a is an operator on meromorphic families in a in which pD p is re- 
placed by a and multiplication by p translates the imaginary part. Since, as 
remarked above, the coefficient of D 2 in E is of the form 0(p 2 ) + O(pv), i.e. 
is a sum of terms having better decay either in the sense of v or p than the 
rest of the operator, we have the following result on the mapping properties 
of R a : 

Lemma 9.1. For each p,u,k,£,s with Imp > — v the operator family R a 
enjoys the following mapping properties: 

(1) R a enlarges the domain of holomorphy: 
(9.1) 

R a :(a - p)- e (H(C u ) n {a)~ k L°° L 2 (R; &\n*S + ))) 
-K(T -{ll- i))- £ {U(C u+l ) n (a)- k+2 L°°L 2 (R; I^(N* S + ))) 

+ {a-{p- 2i))- e (n(C u+2 ) n (a)- k+2 L°°L 2 (R; I^~ 2 )(N* S+))) 

(2) If f a vanishes near C_ for Ima > C, then R a f a also vanishes near 
C- for Ima > C - 1. 

(3) If (p is holomorphic in Ima > C , smooth on X, and rapidly vanishing 
for I Recr| — > 00 in any strip, then 

Ra (V - a )- k - l v ± ^(\og(v ± i0)) fc <; 

k 

= {a- (a - i))-*- 1 ^- 1 ± *0))^, x 

k 

+ {a- (a - Sh))-*- 1 ^" 8 £(log(t; ± *0))^, 2, 

where <j)j t i enjoy the same same properties and are holomorphic on 
Ima > —C — 1. 
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Remark 9.2. Note that in the first term the Sobolev order has decreased 
by 1 arising from the action of the 0(p)vD 2 term in □ (rather than by 2 
as would be the effect of a 0(p)D 2 term). In the second term, we see the 
action of 0(p 2 )D 2 terms, which give a family holomorphic in an even larger 
strip, at the cost of further worsening of Sobolev regularity. We also lose 
at high frequency owing to the (pD p ) 2 error term in the rescaled □, which 
Mellin transforms to an 0(a 2 ). m 

Now we Mellin transform the identity Lu = g, splitting up L = N(L) + E 
to obtain 

Pa(u a ) =9a- R<rU a , 

where, as above, P a = N(L). By compact support in M°, we have for all 
C,s 

~g a G (H(C c )n(a)- co L 00 L 2 (R;I ( - s \N*S + ))). 
As discussed above we know that w lies in some H^ '" 1 (M). Thus by 



Lemma 2.2 
(9.2) 



a* en{c, )n(a) 



max(0, 



-so) L oo L 2 ( 



;H S °) 



where 

ft = 7 - (n - 2)/2. 
(Note that we may nee d to reduce so to guarantee that so + 7 < 1/2 for 

below.) Further, as u vanishes near C_, u a has the 



4.6 



the use of Corollary 
same property. Thus, in the notation of Section [5j the inhomogeneous part 
is in 3^ ftr_1 and u a G X SitT . Here we may choose Sft r to be is constant on the 
singular support of u a as u a is smooth near C_; indeed we may take it to 
be constant except in a small neighborhood of C_ where we have assumed 
w trivial. We take Sf tr to equal sq in the remainder of the space. 



Since we are assuming a non-trapping metric, Corollary 4.6 further tells 
us that w is conormal with respect to p = 



u a G H(C C0 ) n (a)-°°L c 



= 0. By Lemma gj} 

>L 2 (R;I(- 00 \N*S + ))). 



Thus, by interpolation with (9.2), 

u G H(C, ) n (a}- co L co L 2 (R;I (so -°\N*S + ))). 



16 



We further note that a sharper result is true in which the RHS of (9.1 1 is replaced 



by 



f i~i))-\U(C v+1 )n{a)- k+1 L ao L 2 (BL-,I (a - 1) (N*S + ))) 



-k+2 jjx j2 



(n-2i))- l (H(C„+ 2 )n(a}- k L 00 L 2 (R;I <s - 2) {N*S + ))) 
(li-2i))- e (H(C l/+ 2)n{a)- k+1 L 0O L 2 (R;I (s - 1) (N*S + ))) 
(ti-2i))- e (U(C„ +2 )n{a)- k+2 L oc L 2 (R;I M (N*S + ))); 
however, this refinement will not be needed for our argument. 
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By Lemma |9.1[ then, 

R a u a en(C, 0+1 ) n (a)- oo L oo L 2 (R;I^- 1 - \N*S + ))) 

+ n(c, 0+2 ) n (ct)- 00 ^ 00 ^ 2 ^;/^- 2 - ))), 

hence P a u lies in this space as well. 

Because P a u a is now known to be holomorphic in a larger strip, we can 
now invert P a to obtain meromorphy of u a on a larger half-plane - by 



Propositions 5.1 and 5.2 P a is Fredholm as a map 

and -P" 1 has finitely many poles in any horizontal strip Imz G [a, b], and 
satisfies polynomial growth estimates as \Hez\ — > oo. Thus u a is the sum of 
two terms. One is meromorphic in C f0 +i with values in (a)~°° L 2 (Jst, H s °) 
with (finitely many) poles in this strip, arising from the poles of P~ , while 
the other is meromorphic in C ?0 +2 with values in (<7)~°°L 2 (IR, .£P 0-1 ), again 
with (finitely many) poles in the strip. Here (and below) we are ignoring the 
distinction between X Sfti and H s as u a is trivial by hypothesis on the set 
where the regularity in the variable-order Sobolev space differs from H s . 

Since u a enjoys conormal regularity with respect to S+, and (as will we 
will see below) so do the pole terms, we may interpolate the conormal esti- 
mates with the estimates in H s ° (resp. £P 0_1 ) to obtain^] 



+ E ( a - a - 



in ■ 



where 



aj E H(C (o+1 ) n (ay^L^L 2 ^5/(^^+1/2-0)^ _ 
As the aj are given by the polar parts of P^ 1 at values aj lying in a strip in C, 



the coefficients of the polar part of the sum are described by Proposition 8.1 
-l 



3 



K=0 1=0 

(9.3) 



+ 0({a - aj 

with aj K £± = ^^-^)j >m ._^ K ^^± and the 4>j, a ,± are smooth on X. (If — iaj is 



a positive integer, then the term is of the form in Proposition 8.1 ) We may 
further arrange that aj K £± are holomorphic and rapidly decaying in strips. 



1 ^Note the improvement in the Sobolev orders: applying P a 1 we win back the deriva- 
tive we lost from applying R a . 
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Now we iterate this argument: by Lemma |9.1[ 



(9.4) 



R a u e n(c qo+2 ) n {a}- 


" 00 L 00 L 2 (R; l( S()_1_ 


-0))) 


+ n(c, 0+3 )n(a)- 


' 00 L 00 L 2 (M; j( So_2 " 


-0))) 




' 00 L 00 L 2 (M; /( s °~ 3 ~ 


-0))) 


+ E (* 


_( <Ti - l ))- m ^ 




Imo-j>-?o-l 






+ E (* 







where 



and the coefficient of the polar parts of the sums a'j , a" have the form 

m,-l K 



k=0 £=0 



and 



mi— 1 



k=0 £=0 



+ ^->g(<;-*0))^ 



+ 0((a - fa - 2t)) m 

and the bj^-t- and are also smooth (though their values have chang 
Now again inverting P a and employing Proposition 8.1 yields 

u eft(C Co+2 ) n (a^L^L 2 ^; j( s °-°))) 

+ ^(C, 0+3 ) n (a)- 00 ^ 00 ^^; /( s o-i-0))) 

+ H(C Co+4 ) n (a)- 00 ^ 00 ^^; /( s o-2-0))) 

+ E 

+ E (^-(^-or^Sjj 

+ E (a-fa-*))"*^ 

Imffj>- co — 2 
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where again the coefficients of the poles have expansions as in equation (9.3 ) 
(although the expansion for a,j2 begins at v~ l(7j ~ 1 )] here we may have rhj > 
rrij in case there are integer coincidences among the poles of P^ 1 , i.e. if Oj 
and (Tj — i or Oj — 2i are both poles. 
Iterating, we obtain after N such steps: 

u G n(C, 0+N ) n (ay^L^L^R; I^" ))) 
+ ■■■+ n(C f0+2N ) n (ay^L^L^R; I^-N-o)^ 

+ E (<r-<rj)~ mj *J 

Im<Tj>— so— 27V 
N 

+ E E - to - ^r^jr, 

l=\ Im<T j >-c -2Af 

here again fhjt may exceed rrij in case of integer coincidences among poles 
of P~ l . Moreover, while a,j is described by (9.3) above, we also have 

rhj— l i-l P(j,l,n,k) 

&j£ = E ( CT ~ to ~ E E ( v +iO J ~ k ( l °g( v + zO)) p a ifc fcp+ 

k=0 k=0 p=0 

+ vZlo J ~ k (log(v - iO)) p a jeKkp _ 
We now use that inverse Mellin transform of (a — 0"o)~ m is 

„m 

-p-o (1 )m-l 

[m — \y. 

to conclude that, under inverse Mellin transform with a contour deformation 
to the line R — z(?o + N) the poles in the penultimate sum yield the residues: 

( 9 - 5 ) E E E (rr^Hw/^ ( bg Pr— 1 Hoe 
j k=o e=o ^ 3 '' 



i.e. the main terms in our asymptotic expansion q Rearranging this sum 
then shows that it is in fact equal to 



m.i—1 , -, 
,-p+l 



j k=0 p ' 

i.e., the only logarithmic terms in this sum are powers of logu — log p. 



Strictly speaking, this is an expansion in powers of (v ± lO) (or H^v^" 7 if — la is 
a nonnegative integer) rather than in powers of v; however, we are primarily concerned 
with asymptotics in the regime v/p ^> 0. 
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The terms in the last sum become 

N 

EE E a jlKaP ^v-^-^\\ogp\ K \\o g vr, 

j 1=1 K+a<m,j£ 

(which will be lower order following the radiation field blowup, hence we 
do not pursue the question of their structure in much detail). The other, 
"remainder," terms in the sum lie in 

N 

^HQCn+N+j) n ( ff >-°°L°°L 2 CR; i (s °- j - 0) )). 

3=0 



Thus by Lemma 2.2 following inverse Mellin transform our error terms 
become 

Q^0+N+j-0 v s o -l/2-j-0y 

In the radiation field blow-up, these will be 0(p»+iV+ S o-i/2-O s s -i/2-i-0)_ 
Thus, returning to the solution w to D g w = f £ C£°, we find that near 
{p = v = 0}, w has an asymptotic expansion of the form 

w ~ p^r ^ E P™ j v~ ia i (log v - \ogp) K a jK 

j K<rrij 

N 

+p m/2 EE E w^ +e v-^- e+ >^nogv\ a 

j i=\ K+a<rhji 

+ 0{p Co+N+j -°v s °- 1/2 - j - ). 

Note that the only log terms in the p w ii)~ %<J i term occur as powers of log v — 
log p. Because log v — log p = log s in the radiation field blow-up, this implies 
that p~( n ~ 2 ^ 2 w has a restriction to the front face. 

10. THE ASYMPTOTICS OF THE RADIATION FIELD 

The results of the last section yield that solutions of Hw € C C °°(M°) satisfy 
w=p (n-2)/2^2 £ v-^(logv- log p) K a jK 

j K<lllj 

N 

+/"- 2)/2 EE E ^w^ + v^-^ i |iog P riiog«r+^ 

j i=\ K+a<rhji 

with 

N 

w > G y^ /> (n-2)/2+c +Ar+j-O u so-l/2-j-0 L oO j 
3=0 

where aj are the poles of the meromorphic inverse (i-^) -1 , and the co- 
efficients are the corresponding resonance states. Here v~ %a ia 3K (and its 
counterpart in the second sum) is understood to mean a sum of the two 
(v ± z0) -lCTj terms (which we write out fully below). 
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We now introduce the "radiation field" coordinates p,y,s = v/p; note 
that these constitute a coordinate system on the blown up space described 



in Section 3.7, and note that d s is well-defined as a vector field on the fibers 
of ff . In these coordinates, then, homogeneity yields the expansion 

J2 ( 1 °g s ) Q ( a W,+ s +*o J + a' jKai _sZlo J ) 

j a+K<m,j 
N 

+ Y //|log/)| K (logp + logs) a ( 

j £=1 n+a<rhjt 

+ u' 

_ re — 2 

for u = p 2 w. Consequently, restricting terms of the expansion to p = 
yields an expansion 



Y( a i 



-tOj —lOj ~ \ 

ijk,+s +t0 -\-ajk,-S_ t0 -\-Ujk) 

3 

with a remainder term v! . Notice that the presence of logp factors in the 
p° (£ = 0) terms would prevent the restriction of u to the front face of the 
blow-up, but in Section [9] we showed that in fact (at top-order) those terms 
possessing a logarithmic factor cancel. We can now define the radiation field 
as in Section I37T1 



Definition 10.1. If w is a solution of O g w = f, f E C£°(M°), w vanishing 
near C_, we define the (forward) radiation field of w by 



-2 



■fc+M^y) = d s u(0,s,y), u = p 2 w . 

The rest of Theorem |1.1| follows immediately. As identified in Section [7j 
the exponents Oj are the poles of 1Zc+ (o~), i.e., the resonances of the asymp- 
totically hyperbolic problem on the cap C+, while the terms supported at 
S+ do not contribute to the expansion as s — > oo. 

Remark 10.2. While it may seem that the coefficients in the expansion are 
singular at s = 0, this is an artifact of the basis chosen. The b-regularity 
established in Section [4] (see, in particular, Remark 4.7) implies that the 



solution is conormal to the front face of the radiation field blow-up and 
hence the coefficients may be taken to be smooth. 

10.1. Asymptotically Minkowski space. We now consider the special 
case of asymptotically Minkowski space (i.e., "normally very short range" 
perturbations of Minkowski space). Here we are assuming that the metric 



takes the form (3.1) modulo 

0(p)Q 



dp dv duo \ 

P P P J 



with Q a rank- two symmetric tensor with smooth coefficients. Then the 
induced metric on C+ (which is then diffeomorphic to a ball) is the metric 
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on (n — l)-dimensional hyperbolic space; P a is a conjugate of the spectral 
family on hyperbolic space. (See Section 5 of [201) In particular, the relevant 
poles of P^ 1 (i.e., those of L~ + from Section 7| are given by the poles of 

the meromorphic expansion of ^A H n-i — a 2 — ^ n ~^ \ . These poles can 



4 

be calculated explicitly: when n is even (and hence the spatial dimension 
is odd), there are no poles, while if n is odd, the poles are given by Uj = 
— — ij for j £ N. In particular, 7£_|_[u;] has an asymptotic expansion of 
the following form: 



K+[w](s,uj) 




n even 
s~ 2~ J (log s) K dj K n odd. 



(Recall that one differentiates the restriction of u in s to obtain TZ+.) In the 
special case when the metric is in fact exactly Minkowski in a neighborhood 
of C+ in M, we remark that the whole iterative apparatus of Section [9] can 
be dispensed with, in favor of a single application of P" 1 to the Mellin- 
transformed inhomogeneity, with the result that the the log terms in the 
expansion do not appear in that case. 

The stability of P" 1 under perturbations implies that for small "normally 
short range" perturbations of Minkowski space, the radiation field still de- 
cays. In this setting, however, poles of P" 1 that are not poles of L~ , (and 
hence do not affect the decay of the radiation field) may become relevant 
under perturbations. As discussed earlier, such poles must occur at purely 
imaginary negative integers and the corresponding states must be supported 
exactly at S+. Such a state occurs even in 4-dimensional Minkowski space 
at a = —i. Under small "normally short range" perturbations, then, the 
first pole occurs close to a = — i and so we conclude that the radiation field 
is 0(s~ 2+e ) as s — > oo. 



Appendix A. Variable order Sobolev spaces 

First recall that (uniform) symbols a S S™ s on W 1 x IR n of type (p, 5) of 
order r£l are C°° functions on M™ x such that 

\D a z D^a\ < C(C) r+<5|a| ~ p|/31 . 

For various applications, the natural type is p = 1—5, 5 £ [0, 1/2), with 5 = 
corresponding to the standard symbol class. We assume these restrictions 
from now on; for us the relevant regime will be 5 > but arbitrarily small. 
Note that S~°° = f] r S[_ s s is independent of 5. There is a symbol calculus 

within this class S^_ ss , which works modulo S±Zs~s 2S ', the principal symbol 
of the composition of two operators is the product of the two principal sym- 
bols in this sense. Further, one has the full symbol expansion of the compo- 
sition modulo namely if (Au)(z) = (27r)~ n J e l ^ z ^ z ^ ( 'a(z, ()u(z') dz' is 
the left quantization of a £ S[_ s s , and B is the left quantization of b £ S^_ 5 $ 



56 



DEAN BASKIN, ANDRAS VASY, AND JARED WUNSCH 



then AB is the left quantization of 

M 

a 

As usual these can be transferred to manifolds by local coordinates, allowing 
the addition of globally C°° kernels as well. 

We can now turn to variable order operators. Suppose that s is a real- 
valued function on S*R n = K" x S™" 1 = R n x (M n \ {0})/M + , which we 
assume is constant outside a compact set since we are interested only in 
transferring the result to manifolds via local coordinates - one could assume 
instead uniform bounds on derivatives on IR n x S n . On IR n x IR n , we say 
that a is a (variable order) symbol of order s, written a G Sf_ ss , 5 G (0, 1/2) 
if 

(A.l) a = (C) s a , a G S?_ w (r*X). 

So Sl_ ss C S^_ ss with so = sups. Thus, one can quantize these symbols, 
with the result, Vl^^ being a subset of ^^.gg- One calls the equivalence 
class of a in Sf_ ss /S^Zs"s 2S the principal symbol of the left quantization A 
of a. We could of course just as well used another choice of quantization 
such as right- or Weyl-quantization. Note, though, that the condition 5 > 
is crucial for making the different choices of quantizations equivalent since 
the right reduction formula is 

(-z)H 

and the derivatives falling on the exponent of (C) give logarithmic terms, 
which do not have the full S\$ type gain. 

The full asymptotic expansion for composition shows that if s, s' are real 
valued functions on then 

A G *U,6, B G *C S ,6 => AB G 

and modulo it is given by a quantization of the product of the 

principal symbols; again 5 > is important. The commutator [A, B] is 
then in ^if^ 2 ' 5 , and its principal symbol (modulo <S 1 +;f 5 ~ 2+4<5 ) is \{a, b}, 
where {.,.} is the Poisson bracket, and a, b are the respective principal 
symbols. Defining a G S?g S to be elliptic if there exists c, R > such 



E 



— D ° D < a ' 



that \a\ > c(() s for (Q > R, i.e. if ao is elliptic in (A.l) in the analo- 
gous standard sense, the (microlocal) elliptic parametrix construction works, 
i.e. if A G ^{_g g has elliptic principal symbol then there is G G ^7— 8 5 
such that GA — Id, AG — Id G x I / ~ 00 . We can transfer these operators 
to manifolds X via localization and adding C°° Schwartz kernels to the 
space; here we may assume that X is compact. In this manner, for s a 
real- valued function on S*X = (T*X \ o)/IR + with sq = sups, we define 
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*!-<S<sC^) c ^T-ss( X )- Tne Principal symbol of A E ^l_ ss (X) is a well- 
defined element of ' S^_ StS (T*X)/S{Z^ 2S (T*X). 

We can now define Sobolev spaces: fix A E ^ S (X) elliptic, s± = inf s. We 
write 

-ff s = {U E F S1 : AU E L 2 }, ||17|&. = \\U\\jj S1 + \\AU\\ 2 L2 ; 

this is a Hilbert space and all the standard mapping properties of ps.d.o's 
apply. Different elliptic choices A,B£ 5 defining H s give the same 
space, since if G is a parametrix for A, then BU = BGAU + EU, where 
E E so BG E Af7 E L 2 shows BC7 E L 2 by the standard 

L 2 -boundedness of ^i_gg, and also shows the equivalence of the norms. 
Further, if s, s' E C°°(S*X) and B is order s then 

B :H S ' ^ H s '- S 

is continuous; taking A s elliptic of order s, then this is equivalent to 

A s'-s BA -s . L 2^ L 2 

bounded, but the left hand side is in so this is again the standard 

I? boundedness. 

Since the elliptic parametrix construction works, elliptic estimates hold 
without conditions on s in this setting. In our considerations, near the radial 
sets s will be taken constant, so the previous results apply microlocally there. 
However, one needs new real principal type estimates; these hold if s is non- 
increasing along the direction of the H p -flow in which we want to propagate 
the estimates. 

Proposition A.l. Suppose that P E ty m (X) has real principal symbol. 
Suppose that s E C°°(S*X) is non-increasing along H p on a neighborhood O 
ofq£ S*X. Let B,G,R£ with the property that WF'(B) C Ell(G) and 
such that if a E WF'(-B) n S then the backward (null-)bicharacteristic of p 
from a reaches Ell(i?) while remaining in EU(G) n O. Then for all N there 
is C > such that 

\\BU\\ H s < C(\\GPU\\ H3 - m +i + \\RU\\ H s + \\U\\ h -n). 

A similar result holds if s is non-decreasing along H p and 'backward' is 
replaced by 'forward'. 

Related results appear in [T7j , but there the weights arise from the base 
space X, and logarithmic weights are used as well, which would require some 
definiteness of the derivative of s along H p that we do not have here. 

Proof. As the result states nothing about radial points, one may assume 
that Hp is non-radial on O. This then reduces to a microlocal result, namely 
that there is a neighborhood of a point q in which the analogous property 
holds. This can be proved by a positive commutator estimate as in 
Let |£| be a positive homogeneous degree 1 elliptic function on T*X; since 
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we are working microlocally, we may take |£| to be the function |£| in local 
coordinates. With H Pim = |£| -rn+1 Hp denoting the rescaled Hamilton vector 
field, which is homogeneous of degree zero, thus a vector field on S*X, one 
can introduce local coordinates q±, ... , q 2n -\ on S*X centered at a such that 
H p , m = one writes q' = (q 2 , . . . , q 2n -i)- Then one fixes t 2 < h < < t 
and a neighborhood U of in IR 2 ™ -2 such that [t 2 , to] qi X U q > C O and such 

that one has a priori regularity near [t2,ti] gi x U q ', i.e. R in the notation of 
the proposition is elliptic there. For r G [0, 1] (the regularization parameter) 
one considers 

a r = |C| s - (m - 1)/2 x(9i)0(9')^(|CI), 
where G C£°(IR 2n_2 ) is supported in f7, 

V; r (t) = (l+rt)- 1 , 

and 

x(Qi) = Xo(?i)xi(9i)> 
with xo(*) = e~ F ^ t °- t \ t < t , Xo(*) = for t > t and xi(*) = 1 near 
[ii,oo), near (— oo,<2]; here F > will be taken sufficiently large. Taking 
5 G (0, 1/2) arbitrary (i.e. 5 can be very small), ip r reduces the order of a r 
for r > 0, so a r G S^Z^s 1 for r > 0, and for r G [0, 1], a r is uniformly 

bounded in 6**_^ 1 ^ 2 , converging to ao in SlZ^g 1 ^ 2+e for e > 0. Then as 
H p , m <7i = 1 and = n/> 2 , 

H P a* = m 2s 4>(qTM\Z\) 2 Xi(qi)xo(q') 

x (xo(9i)Xi(9i) + Xo(9i)Xi(9i) 

+ (a - (m - l)/2 + r|e|V'r)kr 1 (H p , m |e|)xo(9i)xi('?i) 

+ (Iog|^|)(H P)ro s)xo(9i)xi((?i))- 

Now, Xo < 0, giving rise to the main 'good' term, while the x'i term is 
supported in (t 2 , h) qi x U q >, where we have a priori regularity and estimates. 
Further, by making F large, taking into account that r|£|?/v is bounded, we 
can dominate the |£| -1 H PiJn |£| term since xo can be bounded by a small 
multiple of x'o f° r F > large, and H Pjm s < 0, i.e. has the same sign as 
the x'o term. The imaginary (or skew-adjoint in the non-scalar setting) part 
of the subprincipal symbol also gives a contribution that can be dealt with 
as the |£| -1 H Pim |£| term. Thus, taking A r to have principal symbol a r and 
family wave front set WF'({^4 r }) = esssuppa (for instance a quantization of 
a r ), B r have principal symbol 

b r = \C\ s <P(q')xi(qi)^Xo(q')Xo(q')MO, 

and similar WF' one obtains an estimate of the desired kind, and by esti- 
mating the x'i term (which is the only one having the wrong sign) by the 
R term, first by obtaining an estimate for r > and then letting r — > to 
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obtain the result of the desired form. Corresponding to the symbol class, 
this can give 1/2 — 5 order of improvement (i.e. allows — N = s — 1/2 + 5) 
for all 5 > 0; iterating gives the stated result. □ 
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